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PREFACE 


This  investigation  was  sponsored  by  the  Eustis  Directorate  of  the 
U.  3.  Army  Air  Mobility  Research  and  Development  Laboratory  under  Contract 
DAAJ0P-YU-C-OO27 , Project  1F262209AK76.  The  development  of  the  components 
of  the  Rotorcraft  Wake  Analysis  was  conducted  at  the  United  Technologies 
Research  Center  (UTRC ) over  a period  of  more  than  10  years.  The  Rotorcraft 
Wake  Analysis  was  furnished  as  a result  of  the  on-going  UTRC  research  pro- 
gram. In  accordance  with  the  immediate  objective  of  the  contract,  the 
analysis  was  furnished  for  use  with  the  Rotorcraft  Flight  Simulation  Anal- 
ysis (C-Sl  Program)  available  to  the  Government.  The  correlation,  demon- 
stration, and  documentation  of  the  Rotorcraft  Wake  Analysis,  as  reported 
herein,  were  conducted  under  the  contract. 

The  technical  representatives  of  the  Contracting  Officer  were  Mr. 

Edward  Austin  and  Lt . Lane  Ellington  of  the  Technology  Applications  Divi- 
sion, Eustis  Directorate.  Their  assistance  in  the  demonstration  of  the  com- 
puter program  for  the  Rotorcraft  Wake  Analysis  on  the  Eustis  Directorate 
computer  is  gratefully  acknowledged.  Contract  management  at  the  United 
Technologies  Research  Center  was  conducted  by  Mr.  Anton  J.  Landgrebe 
(Supervisor,  Rotary  Wing  Technology)  and  Mr.  Marvin  C.  Cheney,  Jr.  (Chief, 
Aerodynamics) . 
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INTRODUCTION 


The  task  of  predicting  the  flow  field  induced  by  a helicopter  rotor 
continues  to  be  of  primary  importance  for  providing  and  evaluating  improved 
rotor  designs.  The  continuing  development  of  more  accurate  aerodynamic 
methods  for  predicting  the  temporal  and  spatial  variation  of  velocities  at 
the  rotor  blades  and  in  the  rotor  flow  field  is  resulting  in  significant 
advances  in  rotor  aerodynamic  methodology  important  to  nearly  every  area  of 
rotary -wing  design.  These  areas  include  performance,  vibrations,  structural 
reliability,  and  acoustics.  For  example,  in  addition  to  performance  improve- 
ment, improved  capability  to  predict  the  detailed  flow  characteristics 
results  in  the  prediction  of  more  accurate  blade  loads  and  thus  a more 
rational  design  of  the  blade  structure.  This,  in  turn,  can  lead  to  reduced 
aircraft  development  and  maintenance  costs  and  improved  performance  because 
of  the  resulting  mere  optimized  design.  Also,  the  capability  to  predict 
rotor- induced  flow  velocities  away  from  the  rotor  is  impoHant  for  the  cal- 
culation of  aerodynamic  interference  effects  at  aircraft  components  such  as 
the  fuselage,  wing,  and  tail  surfaces. 

In  order  to  design  and  evaluate  new  rotor  systems,  it  is  recognized  that 
a reasonably  accurate  method  for  predicting  rotor-induced  velocities  is 
required.  It  should  be  capable  of  predicting  velocities  for  all  existing 
and  contemplated  helicopters  operating  over  the  entire  speed  range.  It 
should  be  cost-effective  and  efficient  to  operate.  It  should  represent  the 
current  state  of  the  art,  be  reasonably  accurate,  but  not  be  overly  complex. 
In  addition,  it  should  have  the  growth  potential  to  incorporate  new  tech- 
nology advances  as  they  are  developed  or  become  cost-effective. 

The  major  technological  advance  in  recent  years  to  rotor- induced  velo- 
city methodology  has  been  the  improvement  in  the  mathematical  modeling  of 
the  rotor  wake.  This  has  become  possible  with  the  advent  of  high-speed  digi- 
tal computers.  Many  efforts  have  been  directed  toward  providing  this  meth- 
odology as  shown  in  the  survey  presented  in  Reference  1.  As  a major  part  of 
a long-term  effort  to  advance  the  technology  associated  with  the  aerodynamic 
and  aeroelastic  behavior  of  helicopter  rotor  blades,  the  United  Technologies 
Research  Center,  UTRC  (formerly,United  Aircraft  Research  Laboratories) , has  for 
the  last  decade,  been  conducting  studies  on  the  influence  of  the  rotor  and 
wing  wakes  on  rotor  inflow  and  dependent  rotorcraft  characteristics.  Seme 
of  the  research  programs  conducted  in  this  area  are  discussed  in  References 
1 through  13.  The  scope  of  this  effort  has  included  both  analytical  and 
experimental  programs,  distorted  and  undistorted  wake  analyses,  single-  and 
multi-rotor  configurations,  conventional  and  compound  helicopters,  and  hover 
and  forward  flight  conditions.  Methods  developed  are  currently  in  use  at 
several  organizations  within  the  Government  and  industry. 


With  the  above  requirements  and  capabilities  in  mind,  an  effort  was 
undertaken  to  consolidate  and  expand  several  of  the  available  UTRC  analyses 
to  form  a single  comprehensive  helicopter  rotor  and  wing  wake  analysis  to 
predict  induced  velocities  for  use  with  rotor  performance  and  airloads  anal- 
yses, such  as  the  Government's  Rotorcraft  Flight  Simulation  Program  (Computer 
Program  C-8l) . The  resulting  analysis  is  entitled  the  UTRC  Rotorcraft  Pre- 
scribed Wake  Induced  Velocity  Analysis,  and  will  hereafter  be  referred  to 
briefly  as  the  Rotorcraft  Wake  Analysis.  It  will  be  shown  that  the  Rotor- 
craft Wake  Analysis  is  applicable  to  single  and  dual  (tandem,  side-by-side, 
and  coaxial)  rotor  systems,  and  both  conventional  and  compound  (winged) 
helicopter  configurations.  It  incorporates  a cost-effective  prescribed  wake 
option  which  allows  for  a selection  between  undistorted,  experimental  and 
analytical  wake  representations.  A flow  chart  of  the  Rotorcraft  Wake  Anal- 
ysis, which  will  be  described  in  detail  in  the  body  of  the  report,  is  pre- 
sented in  Figure  1. 


The  accuracy  of  the  component  methods  of  the  Rotorcraft  Wake  Analysis 
has  previously  been  evaluated  in  conjunction  with  rotor  performance  and  air- 
loads analyses  by  comparing  the  resulting  analytical  performance,  airloads, 
blade  response  and  bending  moments  with  test  data.  Many  of  these  compari- 
sons are  included  in  the  aforementioned  references.  In  order  to  provide  a 
more  direct  assessment  of  the  Rotorcraft  Wake  Analysis  by  comparing  the  out- 
put induced  velocities  with  measured  test  data  from  many  sources,  an  induced 
velocity  correlation  study  was  performed,  and  the  results  are  presented  here- 
in. The  use  of  the  Rotorcraft  Wake  Analysis  in  conjunction  with  the  Govern- 
ment's Rotorcraft  Flight  Simulation  Program  has  been  demonstrated. 


Thus,  the  objective  of  the  research  program  reported  herein  was  to 
develop,  evaluate,  and  demonstrate  the  Rotorcraft  Wake  Analysis.  More  spe- 
cifically, it  was  to  provide  an  analysis  for  computing  induced  velocities 
both  at  and  away  frai  the  rotor,  based  on  prescribed  undistorted  and  dis- 
torted wake  representations  for  hover  and  steady  forward  flight  conditions, 
and  for  a wide  range  of  rotor  configurations  including  dual  rotors  and  a 
wing.  In  accordance  with  this  objective,  the  following  sections  consist  of 
(l)  a description  of  the  Rotorcraft  Wake  Analysis,  (2)  results  of  the  induced 
velocity  correlation  study,  and  (3)  conclusions  and  recommendations. 
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DESCRIPTION  OF  THE  ROTORCRAFT  WAKE  ANALYSIS 


COMPONENT  ANALYSES 


The  Rotorcraft  Wake  Analysis  was  developed  by  incorporating  and 
extending  the  following  three  UTRC  prescribed  wake  analyses  and  a rotor- 
wing  interference  analysis  in  a self-consistent,  single  computer  program: 


UTRC  Prescribed  Wake  Rotor  Inflow  Analysis  (UTRC  Deck  F389) 

UTRC  Prescribed  Wake  Hover  Performance  Analysis  (UTRC  Deck  F456) 
UTRC  Coaxial  Rotor  Prescribed  Wake  Hover  Analysis  (UTRC  Deck  G4l3) 
UTRC  Rotor/Wing  Interference  Analysis  (UTRC  Deck  F109) 


In  general,  these  methods  are  not  redundant  in  that  they  are  directed 
toward  a specific  flight  mode  or  type  of  configuration.  They  have  been  com- 
bined and  expanded  to  provide  flexibility  with  regard  to  fli$it  mode  and 
configuration  in  one  computer  program.  Although  the  capability  of  the  com- 
posite Rotorcraft  Wake  Analysis  has  been  extended  beyond  the  scope  of  the 
individual  analyses  listed  above,  these  analyses  form  the  foundation  of  the 
new  computer  program.  Thus,  a brief  description  of  each  is  presented  below 
prior  to  a description  of  the  interrelationship  of  the  analyses  and  the  re- 
finements included  in  the  composite  analysis.  More  detailed  descriptions, 
applications,  and  correlations  of  the  individual  methods  can  be  found  in 
References  1 through  12. 


UTRC  Prescribed  Waite  Rotor  Inflow 


The  function  of  the  UTRC  Prescribed  Wake  Rotor  Inflow  Analysis  (UTRC 
Deck  F389)  is  to  compute  the  time-varying  circulation  and  inflow  distri- 
butions along  the  blades  that  are  compatible  with  a prescribed  set  of  blade 
section  operating  conditions  and  a prescribed  wake  geometry.  This  analysis 
is  primarily  used  for  single- rotor  forward-flight  conditions,  although  it 
is  also  applicable,  but  not  cost-effective,  to  time  independent  hover  condi- 
tions. The  fundamentals  of  this  analysis  generally  represent  an  extension 
of  the  classical  lifting  line  approach  used  successfully  for  fixed-wing 
aircraft,  which  was  initially  applied  to  the  rotary-wing  problem  by  Piziali 
and  DuWaldt  in  Reference  14.  That  is,  although  the  rotary-wing  problem  is 
much  more  complex  than  the  fixed-wing  problem,  due  to  multiple  blades  and 
the  complex  wake  geometry  and  unsteady  flows  due  to  blade  rotation,  the 
fundamental  aerodynamic  relations  of  blade  circulation  to  lift  coefficient, 
angle  of  attack,  blade  motions,  control  settings,  induced  velocity,  and  wake 
geometry  remain  valid  for  this  analysis. 
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The  blade  section  operating  conditions  (nonaerodynamic)  are  generally 
prescribed  frcm  blade  motion  and  control  information  obtained  from  a sepa- 
rate rotor  performance  and  blade  response  analysis.  The  wake  geometry  re- 
presentation is  either  generated  internally  in  the  computer  program  (undis- 
torted wake  for  forward  flight)  or  prescribed  by  input  or  a separate  wake 
geometry  program.  A flow  chart  showing  the  relation  uf  the  iiflow  analysis 
to  the  rotor  performance/blade  response  analysis  and  the  wake  geometry 
analysis  is  presented  in  Figure  2.  Descriptions  and  results  pertaining  to 
the  development  and  application  of  the  rotor  inflow  program  are  presented 
in  References  1 through  7- 

Briefly,  the  mathematical  model  in  the  rotor  inflow  program  consists 
of  the  representation  of  each  blade  by  a segmented  lifting  line,  and  the 
helical  wake  of  the  rotor  by  discrete  segnented  vortex  filaments  consisting 
of  trailing  vovticity,  which  result  from  the  spanwise  variation  of  bound 
circulation.  The  circulation  of  the  wake  for  each  blade  changes  with  azi- 
muth position  and  is  periodic  for  each  rotor  revolution.  The  blades  are 
divided  into  a finite  number  of  radial  segments,  and  the  induced  velocity 
at  the  center  of  each  selected  blade  segnent  is  computed  by  summing  the  con- 
tributions of  each  bound  and  trailing  wake  segnent.  The  contribution  of 
each  vortex  segment  is  obtained  through  use  of  the  Biot-Savart  law,  which 
expresses  the  induced  velocity  in  terms  of  the  circulation  strength  of  the 
vortex  segnent  and  its  geometrical  position  relative  to  the  blade  segment 
at  which  the  induced  velocity  is  desired.  The  bound  circulation  distribu- 
tion is  determined  by  relating  the  wake  circulations  to  the  bound  circula- 
tions, expressing  the  wake  induced  velocities  in  terms  of  the  unknown  bound 
vortex  strengths  by  means  of  the  Biot-Savart  law,  and  developing  a set  of 
simultaneous  equations  relating  the  bound  circulation  and  local  blade  angle 
of  attack  at  each  blade  segment.  These  equations  thus  involve  the  known 
flight  condition,  wake  geometry,  lift- curve  slope,  blade  motion  and  control 
parameters,  and  the  unknown  bound  circulation  values.  The  circulation  equa- 
tions represent  a matrix  of  terns,  the  number  of  which  (number  of  equations 
and  terms  in  each)  equals  the  square  of  the  number  of  points  on  the  rotor 
disk  (number  of  azimuth  positions  times  number  of  radial  stations).  Solu- 
tion of  these  equations  yields  the  desired  bound  circulation  values,  which, 
when  combined  with  the  appropriate  geometrical  relations  in  the  Biot-Savart 
law,  produce  the  required  induced  velocities  at  or  away  from  the  rotor 
blades. 

A simplified  flow  diagram  showing  the  sequence  of  operations  of  the 
program  is  presented  in  Figure  3.  The  known  section  parameters  (airfoil 
noninduced  velocity  components,  lift-curve  slope,  and  geometric  angle  of 
attack)  are  computed.  The  blade-wake  geometry  is  formulated.  The  wake 
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influence  coefficients  (geometric  coefficients)  at  the  blades  are  calculated 
using  the  Biot-Savart  law.  Numerical  techniques  are  used,  to  solve  the  cir- 
culation matrix.  Finally,  the  blade  inflow  and  induced  velocities  off  the 
blades  are  computed. 

More  specifically,  the  rotor  inflow  program  is  based  on  the  following 
formulation  and  assumptions: 

1.  Each  blade  is  represented  by  a lifting  line  (bound  vortex)  divided 
into  a finite  number  of  segnents  (blade  segnents),  each  having  a different 
circulation  strength  (see  Figure  4).  The  aerodynamic  characteristics  at 
the  centers  of  each  segnent  are  assumed  to  be  representative  of  the  entire 
segnent  over  a finite  azimuth  interval. 

2.  The  wake  is  represented  by  a finite  number  of  vortex  filaments 
trailing  from  the  blade  segnent  boundaries.  Each  filament  is  divided  into 
straight  segments,  the  lengths  of  which  are  determined  by  a specified  wake 
azimuth  interval  which  is  equivalent  to  the  azimuth  interval  of  each  blade 
(see  Figure  4).  The  circulation  strength  of  each  trailing  vortex  segnent 
is  constant  along  its  length,  and  is  equivalent  to  the  difference  between 
the  circulation  values  of  its  adjacent  bound  vortex  segments  when  shed  in 
accordance  with  the  Helmholtz  laws  of  conservation  of  vorticity.  The  cir- 
culation strengths  of  different  vortex  segments  along  a vortex  filament 
vary  in  accordance  with  the  variation  of  the  bound  circulations  with  azimuth 
position. 


3.  Viscous  dissipation  effects  on  the  wake  circulation  strengths 
are  neglected,  in  that  the  circulation  of  a given  wake  segnent  is  constant 
with  time.  However,  the  number  of  wake  revolutions  retained  in  the  analysis 
can  be  limited  to  evaluate  an  abrupt  dissipation  of  the  wake.  Also,  the 
viscous  roll-up  of  the  tip  vortex  can  be  approximated  by  combining  tip  fila- 
ments beyond  a prescribed  roll-up  azimuth  interval.  In  addition,  a vortex 
core  size  is  assigned  to  each  vortex  filament.  Potential  theory  is  assumed 
to  apply  outside  the  vortex  core.  Inside  the  vortex  core  an  option  for  rota- 
tional or  zero  flow  is  provided. 

4.  It  is  assumed  that  the  rotor  is  operating  in  steady-state  flight. 
The  inflow  and  wake  from  each  blade  is  assumed  to  be  periodic  with  blade 
spacing.  That  is,  the  inflow  and  wake  geometry  are  the  same  for  each  blade 
when  at  a specific  azimuth  position. 

5.  Shed  wake  segments  (segnents  normal  to  the  trailing  wake  segnents 
mentioned  above)  arising  from  time,  rather  than  radial  variations  of  bound 
vorticity,  are  not  included  in  the  wake  model.  Although  this  emission 
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technically  violates  a Helmholtz  law,  it  is  believed  that  a more  accurate 
representation  of  the  shed  wake  effects  is  obtained  through  the  use  of  experi- 
mental unsteady  airfoil  characteristics  in  the  analysis.  This  implies  that 
the  primary  effects  of  the  shed  vorticity  are  those  associated  with  the  shed 
wake  immediately  behind  the  blade,  and  thus  can  be  approximated  by  those  of 
a fixed-wing  type  of  wake.  Miller  (Reference  15)  shows  that  this  is  reason- 
able at  least  for  moderate  rotor  advance  ratios.  This  approach  not  only 
permits  a reduction  factor  of  two  in  computer  time  and  prevents  unrealistic 
results  associated  with  a finite  filament  wake  model  bat,  more  importantly, 
permits  nonlinear  unsteady  stall  effects  to  be  included  in  a rational  manner 
in  a blade  response  program  through  the  use  of  unsteady  airfoil  data. 

6.  The  wake  geometry  is  prescribed  from  analytical  or  experimental 
results.  Various  options  to  be  discussed  are  available  for  selecting  the 
representation  of  the  wake  model. 

7.  The  airfoil  at  the  blade  is  assumed  to  be  two-dimensional  (radial 
velocity  components  are  neglected).  For  the  linearized  circulation  solu- 
tions and  associated  calculations,  lift-c’urve  slopes,  stall  angles  of  attack, 
and  angles  of  zero  lift  are  provided  which  vary  with  Mach  number.  These 
values  iuujr  be  based  on  unsteady  aerodynamic  data,  if  available,  and  provided 
directly  from  a blade  response  program.  Below  stall,  the  lift-curve  slope 

is  assumed  to  be  constant.  The  blade  section  circulation  is  limited  to  a 
constant  value  for  each  Mach  number  for  angles  of  attack  above  stall. 

8.  In  the  blade-wake  geometry  calculations  the  blades  are  assumed  to 
be  straight  (rigid  bledas).  However,  flexibility  effects  may  be  included 
in  the  circulation  solution  by  providing  the  necessary  noninduced  velocity 
at  each  segnent  associated  with  flexible  blade  motions  and  controls  from  a 
blade  response  program. 

9.  Tangential  and  radial  induced  velocity  components  are  neglected 
in  the  circulation  solution. 

10.  Small-angle  assumptions  are  included  in  the  '■'irculation  solution. 

11.  The  aerodynamic  interference  effects  of  the  rotor  hub,  fuselage, 
and  tail  rotor  are  neglected.  The  rotor  is  assumed  to  be  operating  out 

of  ground  effect. 

Forward-Flight  Wake  Geometry 


The  rotor  inflow  analysis  requires  that  the  rotor  wake  geometry  be 
specified  in  order  for  circulations  and  induced  velocities  to  be  determined. 
Since  the  wake  geometry  is  prescribed,  the  wake  may  vary  from  an  undistorted 
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wake  model  to  a complex  experimental  or  distorted  analytical  wake  model  with 
tip  vortex  roll-up  and  vortex  core  effects  mathematically  modeled.  Thus, 
there  are  several  alternatives  for  rotor  wake  geometry.  The  least  complex 
is  a classical,  undistorted  wake  which  is  simply  a function  of  the  flight 
condition  and  momentum  inflow  velocity.  The  coordinates  of  this  helical 
wake,  which  is  skewed  in  forward  flight,  are  easily  generated  in  this  pre- 
scribed wake  type  of  analysis  given  the  governing  parameters  which  can  be 
iterated  on  if  desired.  The  generation  of  analytical  distorted  wake  geometries 
require  more  complex  and  operationally  expensive  computer  analyses.  Appropri- 
ate experimental  wake  data  are  certainly  most  desirable  but,  except  for  hover, 
are  not  currently  available  for  most  rotor  configurations  and  forward  flight 
conditions . It  has  been  established  that  the  requirement  for  distortions  from 
the  classical  type  wake  geometry  is  dependent  on  the  rotorcraft  configuration, 
flight  condition,  parameter  of  interest,  ana  accuracy  required.  For  example, 
an  undistorted  wake  geometry  is  often  sufficiently  accurate  for  integrated 
performance  calculations  (thrust,  torque,  etc.)  for  conventional  rotorcraft 
operating  in  high-speed  flight.  However,  for  hover  and  low-speed  condit'  >ns, 
wake  distortions  are  very  significant.  Either  an  undistorted  or  a distorted 
wake  geometry  may  be  used  in  the  analysis.  Distorted  wake  geometry  may  be  from 
analytical  or  experimental  sources. 

In  its  simplest  form,  the  wake  frcm  each  blade  can  be  assumed  to  be  a 
classical  undistcrted  skewed  helical  sheet  of  vorticity  defined  frcm  momen- 
tum considerations.  An  option  is  provided  in  the  computer  program  for  repre- 
senting the  roll-up  of  the  outer  portion  of  the  blade's  wake  into  a concen- 
trated tip  vortex.  To  distinguish  between  the  two  optional  wake  represen- 
tations without  wake  distortions,  the  following  nomenclature  will  be  used 
herein: 


Classical  Wake  Representation  — 


The  analytical  wake  representation  in 
which  wake  distortions  are  neglected, 
and  all  vortex  filaments  are  retained 
for  the  same  number  of  wake  revolutions 
(classical  vortex  sheet  representation 
of  the  wake) . 


Undistorted  Wake  Representation  — The  analytical  wake  representation  in 

which  wake  distortions  are  neglected, 
but  the  roll-up  of  the  outer  vortex 
filaments  into  a concentrated  tip  vor- 
tex is  modeled  by  combining  the  outer 
filaments  into  a single  tip  vortex  fila- 
ment at  a prescribed  azimuthal  distance 
from  the  blade.  A vortex  core  size  is 
prescribed  for  the  tip  vortex  filament. 


♦ 
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The  same  undistorted  inboard  wake 
representation  as  in  the  classical  wake 
representation  is  retained. 


Sample  forward  flight  wake  representations  are  shown  in  Figures  5 through  7. 
In  Figure  5>  a sample  computer  plot  of  the  classical  wake  representation 
for  a two-bladed  rotor  at  a 30-kt  flight  condition  is  presented.  All  vor- 
tex filaments  are  shown.  For  the  "undistorted"  wake  representation,  as 
defined  above,  the  five  outboard  filaments  for  this  condition  would  be  com- 
bined at  15  deg  behind  each  blade  to  simulate  each  rolled-up  tip  vortex. 

In  Figure  6,  the  tip  vortex  filaments  are  presented.  (The  filaments  repre- 
senting the  inboard  vortex  sheets  have  been  omitted  in  this  figure  for 
clarity.  However,  they  are  identical  to  the  inner  filaments  of  the  classi- 
cal wake  representation  shown  in  Figure  5.)  The  coordinates  of  the  undis- 
torted wake  representation  are  obtained  from  the  rotor  advance  ratio,  p,, 
thrust  coefficient,  Gj,,  and  angle  of  attack  relative  to  the  tip  path  plane, 
a.  For  example,  the  top  view  of  the  tip  vortex  filaments  in  Figure  6 
is  obtained  directly  fran  the  helicoidal  path  of  the  blade  tip  as  it  trans- 
lates at  the  velocity  Vcos a and  rotates  at  the  velocity  fiR.  The  side  view 
is  dependent  on  the  wake  skew  angle,  X (angle  between  normal  to  the  rotor 
disk  and  wake  boundary),  which  is  defined  in  the  following  manner: 


tan  X - 


ROTOR  FORWARD  VELOCITY  COMPONENT  IN  PLANE  OF  ROTOR  DISK 


MEAN  ROTOR  INFLOW  VELOCITY  NORMAL  TO  ROTOR  DISK 


X = tan 


Vcos  g 
Vs  in  a - 


To  eliminate  the  necessity  for  prescribing  a classical  or  undistorted 
wake  geometry,  input  provisions  are  made  to  obtain  more  realistic  distorted 
wake  geometries  from  experiment  or  a wake  geometry  analysis.  For  example, 
a distorted  wake  analysis  developed  at  UTRC,  entitled  the  Rotor  Wake  Geometry 
Analysis  (UTRC  Deck  F506),  involves  the  following.  First,  an  undistorted 
wake  model  is  defined  along  with  the  distribution  of  circulation  strengths 
of  the  various  vortex  elements  comprising  the  wake.  The  classical  Biot- 
Savart  law  is  then  applied  to  determine  the  velocities  induced  by  each  vor- 
tex wake  element  at  numerous  points  in  the  wake.  These  distorting  veloci- 
ties are  then  numerically  integrated  over  a small  time  increment  to  obtain 
new  vuke  element  positions.  The  process  of  alternately  computing  new  velo- 
cities and  positions  is  continued  until  a converged,  periodic  distorted 
wak*  geometry  is  attained.  Further  details  of  the  procedures  used  to  com- 
pute *uke  geometries  and  results  are  given  in  References  1 through  9* 


Sample  computer  plots  of  the  distorted  tip  vortex  geometry  are  presented  in 
Figures  7 and  8.  In  Figure  7,  the  distorted  tip  vortex  geometry  is  shown 
for  the  same  30-kt  flight  condition  (p  = 0.068)  as  that  of  Figure  6.  In 
Figure  8,  the  side  view  of  the  distorted  tip  filament  of  one  blade  is  com- 
pared with  the  undistorted  tip  filament  for  a higher  speed  condition  (p  = 
O.lt).  In  Figure  9 through  11,  sample  comparisons  frcm  References  2 and  3 
of  predicted  wake  geometries  with  experiment  are  presented  for  low  and 
moderate  speed  conditions. 

If  available,  experimental  tip  vortex  coordinates  may  be  prescribed  in 
the  analysis.  A sample  flow  visualization  photograph  comparing  the  forward 
wake  boundary  of  a ITTRC  model  rotor  indicated  by  tip  vortex  cross  sections 
with  the  classical  wake  boundary  is  presented  in  Figure  12. 

Circulation  and  Induced  Velocity  Solution 

Following  the  formulation  of  t.ie  wake  geometry  representation  in  the 
program,  the  geometric  relations  between  the  wake  and  the  blades  are  calcu- 
lated, and  the  equations  for  the  blade  circulations  and  induced  velocities 
are  applied.  The  Kutta-Joukowski  law  and  the  Biot-Savart  law  are  the  basic 
relations  used  to  obtain  a closed-form  solution  for  these  parameters.  The 
Kutta-Joukowski  law  relates  the  blade  circulation  and  induced  velocity  dis- 
tributions. The  Biot-Savart  law  relates  the  blade  induced  velocity  distri- 
bution to  the  wake  circulations  and  wake  geometry.  The  term  "geometric 
coefficient"  is  used  to  describe  the  influence  coefficients  in  the  latter 
relation,  which  are  functions  only  of  the  wake  geometry.  That  is,  once 
the  circulations  are  computed,  the  velocities  induced  by  the  bound  and 
trailing  vorticity  of  the  rotor  can  be  determined  through  application  of 
the  Biot-Savart  law,  which  can  be  expressed  simply  in  the  following  form: 

v = f (T,  Wake  Geometry) 

A r(r)  (gc) 


Here  the  induced  velocity  at  a point  on  the  blade,  v,  is  a function  of  the 
circulations  and  wake  geometry,  and  is  proportional  to  the  summation  of  the 
products  of  the  circulation  strength,  T,  and  the  geometric  influence  coef- 
ficient, GC,  of  each  element  of  vorticity  in  the  rotor-wake  system.  The 
geometric  coefficient  is  related  only  to  the  geometry  between  a wake  element 
and  the  point  at  which  the  induced  velocity  is  computed.  Equations  for  the 
induced  velocity  components  are  presented  in  Appendix  A.  Several  assumptions, 
mentioned  above,  for  the  airfoil  lift  characteristics,  velocity  components, 
and  inflow  angles  were  included  in  the  establishment  of  the  circulation 
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A flow  diagram  from  Reference  8 showing  the  primary  input  items,  I 

sequence  of  major  operations,  and  output  of  the  computer  program  is  pre-  I 

sented  in  Figure  15.  As  indicated  in  this  figure,  the  program  is  divided  I 

into  three  independent  parts.  The  first  transforms  the  wake  geometry  input  f 

to  wake  coordinates.  The  second  contains  the  computation  of  the  wake  in-  J 

fluence  coefficients  (gea  oric  coefficien  ,s)  at  the  blade  as  defined  by  | 

the  Biot-Savart  law  and  the  numerical  procedures  for  solving  the  circulation  ] 

matrix  and  associated  induced  velocity  distribution.  In  the  third  part,  j 

performance  chars  cteristics  are  computed.  A brief  description  of  the  tech-  i 

niques  used  in  each  of  the  three  parts  is  presented  below.  ! 

i 

i 

Hovering  Waste  Geometry  j 

Generality  regarding  the  specification  of  the  wake  geometry  was  pro- 
vided in  the  computer  program  to  permit  the  evaluation  of  a wide  variety  ! 

of  wake  geometry  models.  This  was  accomplished  by  requiring  only  general 
cylindrical  wake  coordinates  for  the  wake  segment  end  points  in  the  main 
portion  of  the  program.  The  specific  wslce  model  can  then  be  selected  from 
a series  of  separate  wake  geometry  subroutines  representing  various  methods 
for  prescribing  the  wake  shape.  This  facilitates  the  adaptation  of  improved 
wake  models  to  the  program  as  they  become  available.  It  also  permits  the 
assessment  of  classical  wake  models  such  as  the  Goldstein-type  uncontracted 
wake.  Sample  computer  plots  of  a typical  distorted  wake  representation 
(far  wake  instability  neglected)  and  a classical  wake  representation  for  a 
hovering  rotor  are  presented  in  Figure  16. 

For  the  classical  wake  representation,  the  radial  coordinates  of  the 
wake  segment  end  points  comprising  each  trailing  wake  filament  are  everywhere 
equivalent  to  the  blade  radial  coordinate  from  which  they  originated.  The 
surface  boundary  of  the  resulting  wake  model  is  therefore  cylindrical  in 
shape.  The  axial  velocity  of  all  wake  segment  end  points  is  equal  to  a con- 
stant which  is  an  input  item.  The  momentum  induced  velocity  is  normally 
used  for  this  value.  Each  wake  filament  thus  forms  a helix,  and  wake  points 
at  an  equivalent  azimuthal  distance  from  the  rotor  are  also  at  an  equivalent 
axivU.  distance.  A sample  computer  plot  of  a classical  wake  geometry  is 
included  in  Figure  l6.* 


•^Consistent  with  the  forward  flight  wake  discussed  previously,  a distinction 
in  nomenclature  will  be  made  herein  between  undistorted  wake  representations 
with  and  without  tip  vortex  roll-up.  That  is,  the  wake  representation  with- 
out tip  vortex  roll-up  will  be  called  the  "classical  wake  representation," 
and  the  wake  representation  with  tip  vortex  roll-up  will  be  called  the 
"undistorted  wake  representation." 
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The  requirement  for  a method  employing  an  accurate  distorted  wake  model 
for  hovering  conditions  is  established  in  References  8 and  10. It  was  found  that 
the  rapid  contraction  of  the  slipstream  under  a hovering  rotor  places  the 
vortex  system  sufficiently  close  to  the  rotor  blades  that  it  can  cause  sig- 
nificant changes  in  the  radial  distributions  of  induced  velocities  and  a 
loss  in  hover  penormance  relative  to  that  predicted  by  the  classical  wake 
theories . 

For  hovering  conditions,  systematic  model  rotor  wake  geometry  data  have 
been  acquired  experimentally  at  UTRC.  Sample  flow  visualization  photographs 
are  presented  in  Figure  17.  The  data  have  been  generalized  in  Reference  8 
to  facilitate  the  rapid  estimation  of  wake  geometry  for  a range  of  rotor 
designs  and  operating  conditions.  The  generalized  wake  equations  for  the 
hovering  rotor  are  presented  in  Appendix  C.  A special  computer  subroutine 
has  been  prepared  to  model  the  hovering  rotor  wake  in  accordance  with  the 
generalized  wake  representation. 

Rather  than  input  coordinates  for  each  vortex  element  in  the  wake,  the 
wake  equations  presented  in  Appendix  C are  used  to  define  the  wake,  and  the 
individual  coordinates  are  computed  in  the  program.  Use  of  the  wake  equa- 
tions greatly  simplifies  the  input  requirements  while  retaining  sufficient 
accuracy  for  the  computation  of  induced  velocities  at  the  rotor  blades. 

Figure  16  illustrates  how  the  computer  transforms  the  input  wake  con- 
stants to  coordinate  form  and  plots  the  resulting  wake  pattern.  On  the  left 
side  of  this  figure  are  the  top  and  side  views  of  the  computer  representa- 
tion of  a typical  distorted  wake.  For  clarity,  only  the  wake  from  1 blade 
and  12  vortex  filaments  are  shown.  For  this  wake  model,  5 of  the  l6  avail- 
able vortex  filaments  were  used  over  the  outer  8 percent  of  the  blade  to 
represent  the  tip  vortex. 

The  following  is  a description  of  the  major  features  and  assumptions 
of  the  generalized  wake  representation. 

1.  A smoothly  contracting,  stable  wake  is  assumed.  The  axial 
velocity  of  a wake  filament  is  assumed  to  be  constant  within  specified  azi- 
muth regions;  for  example,  the  first  being  between  the  originating  blade  and 
the  following  blade. 

2.  The  wake  filaments  are  divided  into  tip  vortex  filaments  and  vor- 
tex sheet  filaments  for  which  the  geometry  is  defined  independently  by 
separate  sets  of  wake  constants. 
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3.  The  spanwise  division  between  the  vortex  sheet  and  tip  vortex 
region  is  determined  by  the  requirement  that  the  tip  vortex  filaments  have 
the  same  circulation  sense,  and  one  which  is  consistent  with  a negative 
spanwise  derivative  of  the  final  computed  bound  circulation  distribution 
over  the  tip  region  of  the  blade.  This  results  in  the  radial  location  of 
the  peak  bound  circulation  as  the  dividing  point  at  the  blade  between  the 
inboard  vortex  sheet  and  the  tip  vortex  filaments. 

4.  A provision  is  incorporated  in  the  program  for  approximating  the 
roll-up  of  the  tip  vortex  filament  into  a single  filament.  This  is  accom- 
plished by  truncating  the  inner  tip  vortex  filaments  at  an  input  azimuth 
and  assigning  the  prescribed  tip  vortex  geometry  and  combined  circulation 
strength  to  the  remaining  filament  (Figure  16). 

5.  Each  cross  section  of  the  vortex  sheet  at  a specific  wake  azimuth 
is  assumed  to  be  linear  with  radial  position.  The  radial  positions  of  the 
filaments  representing  the  vortex  sheet,  at  a given  azimuth,  are  assumed  to 
be  linearly  proportional  to  the  radial  coordinate  of  the  vortex  sheet 
boundary.  Beyond  an  azimuth  position  corresponding  to  the  azimuth  spacing 
between  blades,  the  vortex  sheet  is  assumed  to  extend  to  the  tip  vortex 
boundary. 

Circulation,  Induced  Velocity,  and  Performance  Solutions 

Following  the  transformation  of  the  analytical  wake  geometry  represen- 
tation from  input  wake  constants  to  wake  coordinates,  the  equations  for  the 
blade  circulations  arid  induced  velocities  are  applied  in  a manner  similar 
to  that  described  previously  for  the  Prescribed  Wake  Rotor  Inflow  Analysis 
used  for  forward-flight  conditions.  The  primary  equations  for  the  Pres- 
scribed  Wake  Hover  Performance  Analysis  are  presented  in  Appendix  D.  Com- 
parison with  the  forward -flight  equations  of  Appendix  B reveals  that  the 
equations  are  similar  but  simplified  due  to  the  absence  of  azimuth  dependency. 

Following  the  solution  for  the  blade- induced  velocity  distribution, 
the  rotor  hovering  performance  is  determined  based  on  conventional  strip 
theory  and  two-dimensional  airfoil  data. 

In  addition  to  improving  the  hover  performance  predictions  of  a wide 
range  of  full-scale  helicopter  rotors  (Reference  l6),  this  analysis  has 
been  modified  and  successfully  applied  to  predict  propeller  performance 
based  on  generalized  propeller  wake  data  (Reference  17)  • 
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UTRC  Coaxial  Rotor  Prescribed  Wake  Hover 


The  UTRC  Coaxial  Rotor  Prescribed  Wake  Hover  Analysis  (UTRC  Deck  C&13) 
was  developed  to  compute  the  blade  inflow,  airloading  and  rotor  performance 
of  a coaxial  type  rotor  such  as  that  of  the  Sikorsky/Army  Advancing  Blade 
Concept  (ABC)  Demonstrator  Aircraft.  Considering  the  success  of  the  pre- 
scribed wake  method  for  single  rotors  based  on  model  rotor  wake  geometry, 
the  prescribed  wake  approach  was  expanded  to  accommodate  dual  rotors.  The 
program  is  a direct  expansion  of  Deck  F456,  described  previously,  except 
that  provisions  have  been  made  for  an  azimuthal  variation  of  blade  loading 
and  prescribed  wake  geometry.  A comparison  of  the  blade-induced  velocity 
and  airload  distributions  for  a model  rotor  at  one  phase  angle  is  presented 
in  Figure  18.  Significant  differences  are  shown  between  the  results  for 
upper  and  lower  blades  which  were  set  at  the  same  collective  pitch  values. 

A sample  UTRC  flow  visualization  photograph  used  to  establish  the  prescribed 
wake  geometry  is  presented  in  Figure  19- 


UTRC  Rotor/Wing  Interference  Ane 


The  UTRC  Rotor/Wing  Interference  Analysis  (UTRC  Deck  F109)  is  a method 
of  calculating  the  wing-induced  flow  through  a rotor  for  use  in  conjunc- 
tion with  rotor  performance  and  dynamic  response  computer  programs.  The 
analysis  calculates  the  rotor  inflow  resulting  from  an  assumed  wing  vortex 
geometry  consisting  of  a classical  lifting  line  and  a series  of  semi-infinite 
trailing  vortex  filaments.  With  this  method,  it  is  possible  to  account  for 
any  desired  spanwise  wing  circulation  distribution  by  varying  the  vortex 
geometry  and  strength. 


The  mathematical  relation  used  to  relate  the  parameters  indicated  in 
Figure  20  as  influencing  the  rotor  inflow  distribution  induced  by  the  wing 
is  the  classical  Biot-Savart  law,  which  was  used  to  express  the  induced 
contribution  of  each  vortex  element  in  terms  of  the  wing  circulation  and 
rotor/wing  geometry.  The  total  wing- induced  velocity  at  each  point  in  the 
rotor  disk  is  calculated  by  superposition  of  the  individual  contribution 
of  all  the  wing  bound  and  trailing  vortex  elements. 


The  rotor/wing  geometry  is  shown  in  Figure  21.  This  analysis  is  des- 
scribed  in  Reference  12,  where  results  are  presented  for  a model  compound 
helicopter.  Contours  of  the  derivative  of  induced  velocity  at  the  rotor 
with  wing  lift  coefficient  are  shown  in  Figure  22.  A sanqale  comparison 
of  the  predicted  rate  of  change  of  rotor  lift  coefficient  with  wing  lift 
coefficient  from.  Reference  12  is  shown  in  Figure  23. 


. 
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CONSOLIDATION  AND  EXPANSION  OF  THE  ANALYSES 

An  mentioned  earlier,  the  Rotorcraft  Wake  Analysis  was  developed  by 
consolidating  and  expanding  the  previously  described  computer  analyses  in 
a self-consistent  single  computer  program.  A flow  chart  of  the  resulting 
Rotorcraft  Wake  Analysis  is  presented  in  Figure  1.  As  shown,  it  is  mainly 
comprised  of  a rotor  analysis  (the  primary  subprogram),  a wing  analysis, 
and  a dual  rotor  program  control.  Descriptions  of  the  components  of  the 
Rotorcraft  Wake  Analysis  and  its  input-output  are  presented  below.  Des- 
scriptions  of  the  refinements  and  additional  capabilities  provided  in  this 
combined  analysis  are  included  in  the  following  sections. 


Rotor 


Basic  Rotor  Ini 


The  basic  rotor  input  to  the  Rotorcraft  Wake  Analysis,  shown  in  Figure 
1,  consists  of  the  flight  condition,  blade  geometry,  airfoil  data,  blade  con- 
trols-, and  response  input.  The  flight  condition  input  consists  of  the  flight 
velocity,  rotor  tip  speed,  and  shaft  angle.  Since  the  blades  and  wake  repre- 
sentations are  based  on  a finite-element  approach  with  numbers  of  elements 
limited  only  by  the  size  of  the  computer,  distributions  of  blade  design  pro- 
perties can  be  adequately  represented.  For  example,  airfoil  variations, 
twist  distributions,  and  chord  distributions  can  be  modeled.  It  includes 
the  capability  to  accept  blade  geometry  inputs  at  arbitrary  radial  spacing. 
Airfoil  data  are  input  and  an  unsteady  airfoil  data  input  option  is  available. 
From  one  to  eight  blades  can  be  accommodated  with  articulated,  teetering, 
or  hingeless  root  conditions.  Blade  controls  and  response  inputs  (collec- 
tive and  cyclic  pitch,  blade  flapping  and,  if  desired,  bending)  must  be  pro- 
vided from  a rotor  performance/blade  response  program  (e.g.,  the  Government's 
Rotorcraft  Flight  Simulation  Program,  C-8l)  or  test  data. 


Wake  Geometi 


Generality  regarding  the  specification  of  the  rotor  wake  geometry  was 
retained  in  the  computer  program.  This  was  accomplished  by  requiring  that 
the  wake  coordinates  for  the  wake  segnent  end  points  be  stored  for  computa- 
tional purposes  in  a prescribed  format;  thus  no  requirement  was  made  as  to 
how  these  coordinates  are  obtained.  This  facilitates  the  adaptation  of 
improved  wake  models  to  the  program  as  they  became  available.  For  this  pro- 
gram, at  the  present  time,  the  following  options  are  available  for  wake 
geometry: 
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1.  A classical  wake  model  option  is  available  for  hover  and  forward 
flight,  in  which  the  classical  wake  coordinates  are  computed  internally. 


2.  An  undistorted  wake  model  option  (with  tip  vortex  roll-up)  is 
available  internally  for  hover  and  forward  flight,  in  which  the  undistorted 
wake  coordinates  are  computed  internally. 


3.  A generalized  wake  model  is  available  for  hover  conditions,  which 
is  based  on  the  generalized  experimental  wake  equations  of  Reference  8 as 
presented  in  Appendix  C . 


4.  Distorted  wake  geometries  (experimental  or  analytical)  may  be 
prescrib;!  in  coordinate  form  via  tape.  Tip  vortex  coordinates  may  be  input 
directly. 


Options  for  the  number  of  wake  revolutions,  vortex  filament  spacing,  wake 
roll-up  and  . ortex  core  representation  and  size  are  included. 


Rotor  Wake  Influence  Coefficients  and  Circulation  Solution 


Following  the  establishment  of  the  rotor  wake  geometry  (see  Figure  1), 
the  Rotorcraft  Wake  Analysis  calculates  rotor  wake  influence  coefficients 
at  the  rotor  and,  if  applicable,  the  wing  and  other  arbitrary  input 
points  based  on  the  Biot-Savart  law  and  a mathematical  representation  for 
a finite  viscous  vortex  core.  Circulation  strengths  at  the  blade  and  in  the 
wake  are  then  calculated  based  on  the  input  flight  condition,  airfoil  data, 
blade  geometry,  controls,  and  response  (with  or  without  the  effects  of  blade 
flexibility)  as  previously  described. 


Induced  Velocity  Solution 


As  shown  in  the  flow  chart  (Figure  1),  die  final  output  of  the  Rotor- 
craft  Wake  Analysis  is  the  temporal  and  spatial  variation  of  induced  velo- 
cities at  the  rotor,  wing,  and  other  arbitrary  points.  These  induced 
velocities  are  calculated  from  the  blade  and  wake  circulations  and  geometric 
influence  coefficients. 


One  of  the  major  refinements  incorporated  in  the  Rotorcraft  Wake  Analysis 
is  the  added  provision  for  calculating  velocities  induced  by  the  rotor  at 
points  removed  from  the  blades.  Three  components  of  velocity  are  calculated 
at  selected  input  points  at  or  away  from  the  rotor  disk  in  the  tip  path 
plane  and  free-stream  coordinate  systems.  Instantaneous  velocities  are 
calculated  for  the  selected  points  at  times  corresponding  to  rotor  azimuth 
intervals.  For  each  point,  an  average  of  the  ins Cantaneous  velocities  over 
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one  rotor  revolution  is  computed  to  represent  the  time-averaged  induced 
velocity.  Since  steady  conditions  are  assumed,  the  instantaneous  velocities 
at  each  point  are  periodic  and  the  time-averaged  velocity  is  representative 
of  all  time.  The  availability  of  the  above  induced  velocity  provisions  ex- 
tends the  capability  of  the  analysis  to  the  determination  of  rotor  flow 
inter!  f'ence  effects  and  permits  comparison  with  experimental  velocity  data 
measured  in  the  rotor  flow  field. 

Rotor  Performance 

Although  the  primary  objective  of  the  Rotorcraft  Wake  Analysis  is  to 
compute  induced  velocities  for  incorporation  in  a separate  rotor  performance 
analysis,  provisions  for  estimating  rotor  performance  have  been  included. 

The  pe vformance  equations  are  based  on  blade-element  theory  for  rigid  blades 
and  i use  of  tabulated  two-dimensional  airfoil  data.  The  performance  cal- 
culai.  -ons  are  based  on  input  blade  motions  and  controls  (blade  response  is 
not  calculated).  Thus,  for  applications  requiring  rotor  trim,  iteration 
with  a more  sophisticated  rotor  performance/blade  response  program  is  neces- 
sary to  result  in  consistent  predictions  of  rotor  perf oimance , blade  motions, 
controls,  and  inflow  distribution. 

Special  Provisions  for  Hover  Conditions 

For  hover  conditions,  it  is  cost  effective  in  a computer  analysis  to 
take  advantage  of  the  symmetrical  wake  structure  and  the  lack  of  azimuth 
variation  of  the  blade  controls,  airloading,  and  response.  When  the  UTRC 
Prescribed  Wake  Hover  Performance  Analysis  (Deck  F456)  was  incorporated  in 
the  Rotorcraft  Wake  Analysis,  this  feature  was  retained.  Also,  the  wake 
generalization  equations  and  experimental  wake  constants  for  hover  presented 
in  Reference  8 have  been  retained  in  the  computer  program  to  facilitate  the 
input  preparation  for  hover  conditions. 

DUAL  ROTOR  AND  CONTROLLABLE- TWIST  ROTOR  PROVISIONS 

The  Rotorcraft  Wak<;  Analysis  has  the  capability  of  computing  induced 
velocities  for  the  following  dual-rotor  helicopter  configurations:  coaxial 

rotor,  tandem  rotor,  and  side-by-side  rotor.  The  dual- rotor  provision 
resulted  from  an  extension  of  the  Coaxial  Rotor  Prescribed  Wake  Hover 
Analysis  (UTRC  Deck  G4l3)  and  the  Prescribed  Wake  Rotor  inflow  Analysis 
(UTRC  Deck  F389)  for  forward  flight  and  noncoaxial  counterrotating  rotors. 
Provisions  are  included  for  blade  phasing,  interrotor  and  wing  interference, 
different  wake  geometries  for  each  rotor,  and  two  sets  of  input  and  induced 
velocity  output.  No  restrictions  are  placed  on  the  relative  positioning  or 
inclination  of  the  two  rotors.  An  option  is  provided  for  calculating 


iivduced  velocities  at  the  blades  of  each  rotor  and/or  at  arbitrarily  selected 
points  in  a selected  reference  axis  system.  As  indicated  in  Figure  1,  a 
dual-rotor  program  control  is  incorporated  to  control  the  flow  of  the  program 
for  each  rotor. 

The  Rotorcraft  Wake  Analysis  is  also  capable  of  calculating  induced 
velocities  of  controllable-twist  type  rotors  (CTR) . The  time-varying  flap 
control  is  provided  for  through  proper  input  of  chord  and  the  unsteady  air- 
foil lift  coefficient  characteristics  for  the  control  flap  blade  section. 

The  controlled  torsional  deflections  due  to  the  prescribed  flap  controls  are 
accounted  for  in  Figure  1 in  the  blade  controls  and  response  input. 


A modified  version  of  the  17TRC  Rotor  Wing  Interference  Analysis  (UTRC 
Deck  FlO?)  was  incorporated  in  the  Rotorcraft  Wake  Analysis.  As  described 
previously,  this  computer  program  calculates  the  inflow  distribution  at  the 
rotor  due  to  a wing  as  rtpiesented  by  a series  of  semi-infinite  horseshoe 
vortices.  Following  the  wing  circulation  solution,  the  Biot-Savart  law  is 
used  to  determine  the  influence  of  the  horseshoe  vortices  on  the  rotor  in- 
flow at  each  blade  station  and  azimuth  angle. 

As  shown  in  the  Rotorcraft  Wake  Analysis  flow  chart  in  Figure  1,  the 
components  of  the  wing  analysis  subprogram  are  essentially  those  of  the 
Deck  F109  flow  chart  previously  presented  in  Figure  20.  The  analysis  has 
been  modified  and  expanded  to  include  the  following  features: 

1.  The  capability  to  iterate  between  the  rotor  inflow  distribution 
and  wing  inflow  distribution  to  account  for  the  mutual  rotor/wing  wake  inter- 
ference. 

2.  The  capability  to  calculate  the  wing  circulation  and  induced 
velocity  distribution  and  associated  wing  lift  coefficient  in  a self-consis- 
tent iterative  manner.  This  provides  an  alternate  to  the  use  of  Shrenk's 
rule  (as  used  in  Reference  12)  and  also  a means  for  including  the  induced 
velocity  due  to  the  rotor  in  a more  accurate  wing  inflow  iteration. 

3-  The  capability  for  arbitrary  relative  positioning  between  the 
rotor  and  wing. 


4.  The  capability  to  include  inclination  of  the  wing  wake  to  pro- 
vide a first  approximation  for  wing  wake  deflection  as  influenced  by  the 
rotor  wake. 
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The  provision  to  iterate  tetween  the  rotor  and  wing  inflow  distributions 
was  incorporated  by  making  use  of  the  expanded  rotor  analysis  capability  of 
computing  off-rotor  induced  velocities.  Time-averaged  velocities  induced 
by  the  rotor  at  the  wing  stations  are  added  to  the  wing's  self- induced 
velocity  distribution  to  determine  the  wing  circulation  distribution. 
Iteration  between  the  rotor  characteristics,  including  wing  interference 
effects,  and  wing  characteristics,  including  rotor  interference  effects, 
results  in  a self-consistent  solution  for  rotor  and  wing. 


Options  for  prescribing  or  calculating  the  wing  circulation  distribu- 
tions have  been  provided.  If  calculated,  one  approach  is  to  use  the  experi- 
mental wing  downwash  factor  (Reference  12)  to  obtain  the  wing  lift/rotor 
lift  derivative  and  then  to  appropriately  scale  Shrenk's  wing  circulation 
distribution.  This  does  net  account  for  the  spanwise  and  timewise  variation 
of  rotor  interference  at  the  wing.  The  second  approach  is  the  lifting  line 
solution  in  which  spanwise  variations  are  accounted  for. 

Computer  Program  Cost  Efficiency,  Operational  Simplicity,  and 
Growth  Potential 


In  addition  to  accuracy,  a primary  concern  in  the  development  of  the 
several  UTRC  wake  analyses  has  been  cost  efficiency.  Avoidance  of  redun- 
dant and  insignificant  calculations  and  inclusion  of  a fast  iterative  tech- 
nique for  solving  large  matrices  have  resulted  in  significant  reductions 
in  computer  time  requirements.  Although  the  Rotorcraft  Wake  Analysis  is  a 
large  computer  program,  careful  control  of  program  call  sequences  and  overlay 
techniques  has  resulted  in  the  maintenance  of  program  storage  and  computer 
time  efficiency  and  the  allowance  for  the  advantages  of  a single,  composite, 
self-consistent  program. 

Attention  has  been  given  to  operational  simplicity  for  the  computer 
analysis.  Commonality  between  the  input-output  of  the  various  subprograms 
has  been  provided.  Provisions  for  self-consistent  iterations  within  the 
program  have  been  included  along  with  a multiple  case  option. 

Through  the  division  of  the  Rotorcraft  Wake  Analysis  into  component 
subprograms,  growth  potential  has  been  provided,  in  that  any  subprogram  can 
be  updated  or  replaced  as  the  state  of  the  art  advances  or  the  need  arises. 
Also,  the  use  of  a prescribed  wake  in  the  analysis  retains  flexibility  on 
the  type  and  source  of  the  wake  geometry.  As  additional  accurate  and 
generalized  wake  geometry  data  becomes  available,  it  may  be  readily  incor- 
porated in  the  analysis. 
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INDUCED  VELOCITY  CORRELATION  STUDY 


In  order  to  demonstrate  the  accuracy  of  the  Rotorcraft  Wake 
Analysis,  a correlation  study  was  conducted  in  which  the  analysis  was 
applied  to  compare  with  induced  velocity  test  data  from  many  sources. 
The  selection  of  the  test  data,  the  limitations  of  the  test  data,  the 
procedures  used  to  apply  the  analysis  to  the  test  conditions,  and  the 
results  of  the  correlation  study  for  both  hover  and  forward-flight  con- 
ditions are  presented  below. 

SELECTION  OF  INDUCED  VELOCITY  TEST  DATA 


In  the  selection  of  the  induced  velocity  test  data,  it  was  recog- 
nized that,  in  addition  to  being  accurate,  the  data  selected  should 
cover  a wide  range  of  rotorcraft  configurations,  sizes,  flight  conditions 
and  test  point  locations  relative  to  the  rotor(s).  That  is,  it  was 
desired  to  compare  predicted  velocities  with  measured  values  for  points 
near  the  rotor  disk,  and  points  within  and  outside  of  the  rotor  wake 
for  various  combinations  of  flight  speed  and  thrust  level.  Test  data 
that  include  hover  and  low-  and  high-speed  flight  conditions  for  the 
various  types  of  configurations  that  the  analysis  is  applicable  to,  i.e. 
single,  tandem,  coaxial,  side-by-side,  compound  (winged),  and  CTR  type 
rotorcraft,  were  desired. 

To  establish  what  induced  velocity  data  existed  for  selection  in 
this  correlation  study,  a search  of  the  available  literature  on  this 
subject  was  conducted.  Ten  sources  of  data  were  selected  (References 
18  to  27).  Pertinent  information  on  each  are  listed  in  Table  1.  The 
sources  include  data  accumulated  over  a time  span  of  two  decades  (195** 
to  1974)  by  Army,  NASA,  and  industry  organizations. 

The  Rotorcraft  Wake  Analysis  was  applied  to  25  rotor  conditions 
(combinations  of  rotor  configuration  and  test  condition)  selected  from 
the  available  test  data.  Although  the  data  do  have  significant 
limitations  (to  be  described),  the  overall  range  of  the  data  selected 
from  the  references  is  quite  comprehensive  in  that  the  following  are 
included: 

1.  Hover  and  forward -flight  conditions  (g,  = 0 to  0.23) 

2.  Single  and  tandem-rotor  configurations 

3.  Model  and  full-scale  rotors  (R  = 1 to  17.6  ft) 
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4.  Model  hover  test  stand,  outdoor  whirl  tower,  wind  tunnel 
(4  x 6 ft  to  40  x 80  ft)  and  flight  test  data 


5.  Rotors  with  2,  3>  end  7 blades 


6.  Thrust  coefficients  (CT)  from  0.002  to  0.010 


7.  Tip  speeds  (f®)  from  110  fps  to  700  fps 


8.  Varying  blade  designs  (twist,  taper,  airfoil  section,  and  root 
constraint ) 


9.  Local  time-averaged  (mean)  and  local  time-dependent  (instanta- 
neous) velocities 


10.  Three  velocity  components  (some  sources) 


11.  Velocities  measured  by  pressure  (Pitot-static  tubes  and  dynamic 
pressure  disks),  anemometers  (hot-wire  and  hot-film),  and  laser 
velocimeter  techniques 


12.  Velocities  at  points  near  and  far  from  the  rotor,  and  inside 
and  outside  of  the  rotor  wake 


LIMITATIONS  OF  THE  INDUCED  VELOCITY  TEST  DATA 


The  general  limitations  of  the  available  induced  velocity  test  data 


1.  No  one  test  contained  data  which  were  sufficiently  extensive 
to  provide  a single  set  of  accurate,  self-consistent  data  for  ab- 
solute correlation  of  time-averaged  and  instantaneous  velocities 
for  all  desired  locations  over  a wide  range  of  conditions. 


2.  Most  of  the  data  sources  did  not  measure  or  properly  document 
all  of  the  associated  information  required  for  the  Rotorcraft 
Wake  Analysis.  That  is,  in  most  cases,  the  complete  combination 
of  rotor  design  parameters,  pitch  control  angles,  blade  response 
(coning  and  flapping),  and  rotor  performance  (forces  and  moments) 
was  not  known  or  noted  in  the  reference.  Although  sufficient  in- 
formation was  generally  provided  to  estimate  the  missing  informa- 
tion or  to  determine  it  from  a rotor  performance  analysis,  this 
detracted  from  the  absolute  correlation  of  the  Rotorcraft  Wake 
Analysis . 
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3.  In  most  cases,  wake  geometry  was  not  measured  along  with  the 
induced  velocities.  In  no  case  were  blade  airloads  measured.  If 
measured , these  items  and  more  complete  rotor  performance  data 
would  have  been  of  significant  value  for  checking  the  analytical 
wake  representations  and  circulation  solutions  used  to  compute 
induced  velocities. 

Only  one  of  the  sources  of  nover  data  included  quantitative 
instantaneous  induced  velocities.  The  others  were  limited  to 
time-averaged  (mean)  velocities.  The  sources  of  forward- flight 
data  with  meaningful  instantaneous  velocities  reported  on  initial 
demonstration  tests  using  recently  developed  laser  velocimeter 
techniques,  and  they  were  limited  to  selected  points  at  one  or  two 
blade  asimut'n  jositions. 

5.  Except  fur  the  laser  velocimeter  data,  all  available  velocity 
data  were  measured  away  from  the  rotor  blades.  The  measurement 
of  the  instantaneous  induced  velocities  at  the  blades  over  a com- 
; lete  revolution,  which  is  of  extreme  interest  to  the  rotor  de- 
signer, has  yet  to  be  attempted. 

6.  -low  velocity  data  remain  to  be  acquired  for  many  types  of 
rotor  c onfigurati .ns  and  test  conditions.  In  particular,  informa- 
tion is  almost  completely  lacking  for  dual  rotors  (tandem,  coaxial, 
and  side-by-side  rotors),  comj ound  (winged;,  and  controllable 
twist  (CTr.  conf  igurations . Also,  available  velocity  data  for 
forward-flight  conditions  and  varying  combinations  of  blade  designs 
are  very  limited. 
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7.  The  accuracy  of  ti.e  available  velocity  data  is  limited  to  the 
recognised  degree  >if  accuracy  obtainable  with  the  respective  type 
of  equipment  used,  a:. j by  the  sire  of  the  rotor  at.d  the  procedures 
and  facilities  used  by  the  test  conductor,  ror  example.  Pi  tot - 
static  tubes,  iynamic  : ressure  disks,  and  hot-wire  equipment  have 
known  limitation-  regarding  flow  angularity  and  magnitude.  In  most 
sources  of  data,  the  accuracy  of  the  measured  velocities  was  not 
noted.  This  makes  questionable  the  extent  to  which  the  velocities 
predictei  with  the  Rotorcraft  Wake  Analysi.  can  be  confidently 
compared  with  the  test  data  on  an  absolute  basis. 

The  test  data  were  not  always  acquired  under  the  ideal  environ- 
mental conditions  and  rotor  response  behavior  assumed  by  the 
analysis.  The  analysis  is  based  on  the  assumntion  of  steady-state 
conjiti.ns  in  forward  flight  and  the  absence  of  wind  in  hover.  In 
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practice,  particularly  under  realistic  outdoor  conditions,  ideal 
conditions  are  seldom  present,  as  was  attested  to  in  several  data 
sources.  Also,  the  blades  were  assumed,  to  be  manufactured  iden- 
tically and  to  respond  accordingly.  The  sensitivity  of  local  in- 
duced velocities  to  wind  and  blade-to-blade  imperfections  is  re- 
cognized, especially  at  hover  and  low-speed  conditions,  and  for 
some  cases  influence  the  degree  of  velocity  correlation  on  a local 
point  basis. 

Although  the  above  list  of  test  data  limitations  appears  formidable, 
it  is  not  intended  to  imply  that  the  forthcoming  correlation  results 
are  not  useful  for  determining  the  accuracy  of  the  Rotorcraft  Wake 
Analysis.  The  intent  will  be  to  show  that,  although  in  single  cases 
questions  may  arise  regarding  the  absolute  comparison  of  the  predicted 
and  test  results,  in  total  the  accuracy  of  the  analysis  in  its  current 
form  can  be  reasonably  assessed  on  an  initial  basis.  With  the  recent 
advances  in  velocity  measurement  techniques,  it  is  expected  that  an  ex- 
tensive amount  of  accurate  induced  velocity  data  will  become  available 
to  further  evaluate  the  analysis  and  to  provide  direction  for  its 
refinement.  Continued  application  of  the  Rotorcraft  Wake  Analysis  in 
conjunction  with  rotor  performance  and  blade  response  analyses  to  cor- 
relate airloads,  blade  motions,  and  bending  moments  will  also  demon- 
strate the  accuracy  of  the  analysis. 

CORREIATION  STUDY  FOR  HOVERING  CONDITIONS 


Six  of  the  selected  sources  of  data  (References  18  to  23)  contain  informa 
tion  pertaining  to  hovering  test  conditions . Comparisons  of  results  of  the 
Rotorcraft  Wake  Analysis  with  test  data  from  each  are  presented  below. 

Comparison  With  Test  Data  of  Boatwright,  USAAMRJL  TR  72-33 

Induced  velocity  data  for  a full-scale,  two-bladed,  0H-13E  helicop- 
ter rotor  tested  outdoors  on  a 60-ft  rotor  test  tower  are  presented  in 
the  following  reference  (Reference  18); 

Boatwright,  D.W. , Measurements  of  Velocity  Components  in  the 
Wake  of  a rull-Scale  Helicopter  Rotor  in  Hover,  Mississippi  State 
University,  USAAMRDL  Technical  Report  72-33,  Eustis  Directorate, 

U.S.  Army  air  Mobility  Research  and  Development  Laboratory,  Fort 
Eustis,  Va.,  August  1972,  AD  754644. 

This  report  presents  three-component  velocity  measurements  made  below 
the  rotor  with  a splj.t-film  total  vector  anemometer.  Test  data  were 
acquired  by  traversing  the  total  vector  probe  along  an  instrumentation 
boom  which  extended  radially  into  the  rotor  wake. 
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The  boom  was  positioned  at  preselected  distances  below  the  rotor,  and 
instantaneous  velocities  were  measured.  The  term  instantaneous  will 
be  used  herein  in  the  broad  but  conventionally  used  sense  (assuming 
periodicity)  to  indicate  that  the  velocity  is  for  a fixed  blade  azimuth 
position  relative  to  a measurement  point.  This  distinguishes  it  from 
the  mean  or  time-averaged  value  obtained  by  time  averaging  the  velocity 
at  a point  over  a continuous  time  period.* 


Information  on  the  rotor  design  and  test  conditions  are  presented 
in  Table  2.  As  noted  in  the  table,  data  for  three  test  conditions  were 
presented  which  include  independent  variations  of  thrust  coefficient 
(by  varying  collective  pitch)  and  tip  speed.  The  measured  thrust  co- 
efficient values  were  accurately  piedicted  by  using  the  test  collective 
pitch  values  in  the  Rotorcraft  Wake  Analysis.  Predicted  rotor  torque 
levels  for  the  generalized  wake  representation  were  in  the  vicinity  of 
the  high  limit  of  the  range  of  measured  torque  data  presented  in  the 
reference.  Predicted  torque  levels  for  the  undistorted  wake  representa- 
tion were  generally  higher  than  the  measured  values. 

Predicted  induced  velocities  were  oVained  from  six  computer  cases 
of  the  Rotorcraft  Wake  Analysis  comprising  the  combinations  of  the  three 
test  conditions  with  two  types  of  prescribed  wake  representations.  The 
two  wake  representations  used  were  the  undistorted  wake  and  the  general- 
ized wake  based  on  the  model  rotor  wake  data  of  Reference  8 (see  Figure 
17).  The  undistorted  wake  model,  as  defined  earlier,  is  the  same  as 
the  classical  wake  model  except  that  it  combines  the  tip  filaments  of 
the  classical  wake  into  a single  undistorted  tip  filament.  This  wake 
model  was  used  initially  for  this  source  of  data  instead  of  the  classical 
wake  because  it  was  potential. .y  a better  representation  for  correlating 
with  instantaneous  velocities. 


* In  USAAMRDL  Report  72-33,  the  term  instantaneous  is  used  in  the  exact 
sense  to  indicate  a velocity  measurement  taken  at  a point  at  one  instant 
of  time.  The  "mean"  values  in  that  report  are  not  the  continuous  time- 
averaged  values  at  a point,  but  rather  the  instantaneous  velocities  for 
one  rotor  azimuth  position  averaged  for  25  rotor  revolutions,  each  of 
which  would  be  equivalent  if  there  were  no  unsteadiness  of  the  hovering 
rotor  and  wake.  Thus,  they  are  actually  a time-dependent  mean  of  the 
local  instantaneous  velocities  corresponding  to  one  rotor  position. 

Since  wake  unsteadiness  is  not  incorporated  in  the  Rotorcraft  Wake 
Analysis,  instantaneous  velocities  are  calculated  assuming  periodicity, 
which,  although  labeled  instantaneous  herein,  correspond  to  the  "mean" 
values  of  the  reference. 
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TABLE  2.  ROTOR  AND  TEST  INFORMATION  — BOATWRIGHT,  USAAMRDL  72-33 


Rotor  Type 

OH-13E,  single,  full 

Rotor  Design 
Number  of  blades 
Rotor  radius,  ft 
Root  chord,  ft 
Tip  chord,  ft 

Test  Conditions 
Condition  1:  /a  = 0, 
Condition  2:  /a  = 0, 
Condition  3:  /a  = 0, 


scale,  teetering  rotor 


2 

Rotor  solidity 

0.03625 

17.56 

Blade  taper  ratio 

1.38 

1.167 

Blade  twist,  deg 

-4.25  (avg.) 

0.845 

Airfoil  section 

NACA  0015 

flR  = 625  fps,  0 = 6.25  deg,  CT  = 0.0020,  Vq  = 19.8  fps 
HR  = 450  fps , 0 = 10.75  deg,  CT  = 0.0040,  Vq  = 20.1  fps 
fiR  = 450  fps,  0 ^ = 6.25  deg,  CT  = 0.0020,  vq  = 14.2  fps 


Two  uncertainties  were  present  when  the  wake  representations  were 
established.  First,  blade  tip  position  due  to  coning  and  blade  bending 
were  not  indicated  in  Reference  18  and  thus  had  to  be  estimated.  Second, 
the  wake  generalization  of  Reference  8 does  not  include  provision  for 
tapered  blades,  and  thus  the  influence  of  taper  on  the  generalized  wake 
geometry  was  not  included. 

The  positions  of  the  tip  vortices  of  the  generalized  wake  model, 
as  indicated  by  the  cross-section  locations  of  the  tip  filaments,  are 
shown  in  Figure  24.  The  tip  vortex  positions  are  shown  for  the  three 
test  conditions  and  the  following  two  instants  of  time:  (1)  when  the 

reference  blade  is  at  f = 0 deg,  which  indicates  that  the  blade  is 
directly  above  the  measurement  point;  and  (2)  when  the  reference  blade 
is  at  \|f  = 90  deg.  As  shown,  the  generalized  wak'  positions  are  identi- 
cal for  test  conditions  1 and  3 due  to  the  independence  of  hovering 
rotor  wake  geometry  to  tip  speed  variations  at  constant  thrust  co- 
efficient, as  described  in  Reference  8.  Also  indicated  in  Figure  24 
are  the  z/R  elevations  of  the  probe  locations  and  the  approximate  mean 
of  the  measured  radial  excursions  of  the  tip  vortices  at  these  eleva- 
tions. Measured  vertical  and  radial  excursions  of  tip  vortex  points 
as  large  as  5 to  10  percent  of  the  rotor  radius  are  indicated  in  the 
reference,  even  for  near-wake  points.  These  excursions,  due  mainly  to 
ambient  wind  effects,  include  the  generalized  wake  positions  in  the 
near  wake  (at  z/R  = -0.1  and  -0.3  in  Figure  24).  Figure  24  will  be 
useful  in  describing  the  induced  velocity  results. 

Comparisons  of  the  measured  and  theoretical  instantaneous  induced 
velocities  are  presented  in  Figures  25  through  29.  The  induced  veloc- 
ities have  been  nondimens ionalized  by  the  momentum  theory  values  of 
induced  velocity  at  the  rotor,  which  are  19-8,  20.1,  and  14.2  fps  for 
test  conditions  1,  2,  and  3?  respectively.  In  Figure  25,  a comparison 
of  the  radial  distribution  of  the  three  velocity  components  is  pre- 
sented for  test  condition  1 at  the  vertical  location  nearest  to  the 
rotor;  z/R  = -0.1.  The  components  are  based  on  a right-handed  coordi- 
nate system.  Considering  the  plane  of  measurement  points  to  be  at  f = 0 
results  in  the  normal  component,  v , positive  up:  the  radial  component, 
v , positive  outward;  and  the  tangential  component,  v , positive  in  the 
direction  of  rotation.  The  generalized  wake  results  exhibit  good  corre- 
lation with  the  test  data  and  a significant  improvement  over  the  undis- 
torted wake  results.  The  latter  is  most  evident  near  the  wake  boundary, 
where  the  results  are  very  sensitive  to  the  accurate  positioning  of  the 
tip  vortices,  both  radially  (wake  contraction)  and  vertically.  In 
Figure  26,  the  normal  velocity  component  is  compared  at  three  elevations 
below  the  rotor  for  test  condition  1.  Good  agreement  is  indicated 
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between  the  generalized  wake  results  and  the  test  data  at  z/R  = -0.1 
and  -0.3.  At  z/R  - -0.7,  the  correlation  is  poor  near  the  wake 
boundary.  This  is  due  to  the  known  shortcoming  of  the  generalized 
wake  of  Reference  8,  which  was  modeled  therein,  concentrating  on 
near-wake  accuracy  for  predicting  induced  velocities  at  the  rotor.  As 
demonstrated  in  both  the  model  rotor  flow  visualization  results  and  the 
results  of  the  Wake  Geometry  Analysis  in  Reference  8,  the  far  wake  of 
a hovering  rotor  does  not  continue  to  contract  smoothly  as  assumed  in 
the  generalized  wake  representation,  but  rather  expands  in  an  unsteady 
and  unstable  manner.  This  wake  expansion  was  evident  in  Boatwright's 
full-scale  rotor  test  data,  as  shown  in  Figure  24 ; its  effect  is  shown 
in  Figure  dG  at  z/R  = -0.7  by  the  outward  shift  of  the  measured  velocity 
distribution.  An  indication  of  the  capability  of  the  Rotorcraft  Wake 
Analysis  to  more  accurately  account  for  this  if  the  generalized  wake 
were  refined  to  include  the  wake  expansion  is  shown  by  the  outward  shift 
of  the  undistorted  wake  velocities  relative  to  the  generalized  wake 
velocities  due  to  the  complete  lack  of  contraction. 

In  Figure  27,  the  results  for  test  conditions  2 and  3 are  added  to 
the  results  for  test  condition  1 presented  in  Figure  26.  As  stated  in 
Reference  18,  when  the  measured  velocity  components  are  nondimensional- 
ized  by  the  momentum  values  of  induced  velocity,  no  large  differences 
between  the  radial  distribution  of  the  velocity  components  are  generally 
detected.  The  differences  that  do  occur  could  not  be  directly  attri- 
buted to  variations  of  the  test  parameters,  because  variations  due  to 
ambient  wind  effects  produced  wake  geometry  deviations  between  test 
periods  which  masked  any  possible  deviations  due  to  test  parameters. 

This  was  particularly  true  at  the  wake  boundary,  where  unsteady 
variations  of  the  paths  and  transport  velocities  of  the  tip  vortices 
between  repeat  test  conditions  were  found  to  be  large  even  near  the 
rotor.  Also,  a mismatch  of  blade  twist  rate  between  the  two  blades 
(-4.0  and  -4.5  deg)  resulted  in  wake  geometry  variations  between  the 
two  blades. 


The  generally  small  variations  with  test  conditions  of  the  pre- 
dicted velocity  results  from  the  Rotorcraft  Wake  Analysis,  also  shown 
in  Figures  27  and  28,  are  in  agreement  with  the  test  findings.  The  pre- 
dicted results  for  test  conditions  1 and  3 are  identical  (same  general- 
ized wake  was  used),  indicating  that  nondimensionalized  induced  velocity 
is  independent  of  tip  speed.  The  predicted  velocity  variations  due  to 
thrust  coefficient  are  generally  small.  However,  the  existence  of  ex- 
ceptions is  indicated  in  Figure  28  for  the  radial  component  at  z/R  = 
-0.1,  where  the  theory  predicts  a large  local  radial  inflow  at  the  wake 
boundary  for  conditions  1 and  3 which  is  not  predicted  for  test  condi- 
tion 2.  The  presence  of  large  local  velocity  gradients  such  as  this  is 
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explained  by  the  close  proximity  of  a point  to  a tip  vortex.  In 
Figure  24,  it  is  shown  that  when  the  reference  blade  is  at  \|i  = 0 deg, 
a tip  vortex  is  located  just  below  the  z/R  = -0.1  elevation,  which, 
considering  the  counterclockwise  circulation  sense  of  the  tip  vortex, 
produces  the  large  radial  inflow  near  the  0.85  radial  station.  For 
test  condition  2,  the  tip  vortex  positions  corresponding  to  $ = 0 deg 
are  not  near  the  z/R  = -0.1  elevation,  and  thus  a large  velocity 
gradient  is  not  produced.  The  test  data  in  Figure  28  indicate  a large 
radial  velocity  near  the  wake  boundary  for  test  condition  1,  which  is 
consistent  with  the  predicted  results.  However,  unlike  the  predicted 
results,  the  magnitude  is  decreased  for  test  condition  3-  This  is 
probably  attributable  to  the  movement  of  the  tip  vortex  caused  by  am- 
bient wind  effects,  and  illustrates  the  sensitivity  of  local  velocity 
magnitudes  to  tip  vortex  relocations  due  to  environmental  factors. 


The  sensitivity  of  the  local  instantaneous  induced  velocities  to 
tip  vortex  position  is  further  illustrated  in  Figure  29,  in  which  the 
theoretical  velocities  are  compared  with  the  test  data  at  z/R  = -0.1 
for  test  conditions  1 and  2 when  the  reference  blade  is  90  deg  ahead  of 
the  measurement  azimuth  (<li  - 90  deg).  In  Figure  29,  the  measured 
normal  velocities  indicate  the  close  proximity  of  a tip  vortex  near 
x/R  = 0.83  for  both  tes+  conditions.  The  theoretical  results  are  con- 
sistent with  those  of  the  test  for  test  condition  2,  but  not  for  test 
condition  1.  This  is  explained  by  referring  again  to  Figure  24  and 
noting  that  a tip  vortex  in  the  generalized  wake  representation  is  lo- 
cated just  above  the  z/R  = -0.1  elevation  for  test  condition  2 when 
V = 90  deg,  but  not  for  test  condition  1.  The  requirement  for  accurate 
tip  vortex  positioning  is  emphasized  further  in  Figure  29,  where  the 
predicted  radial  velocity  component  for  test  condition  2 is  similar  to 
the  test  results  in  exhibiting  a large  peak  magnitude  at  the  wake 
boundary,  but  is  opposite  in  direction  radially  (outward  rather  than 
inward).  This  is  easily  explained  by  the  generalized  wake  position  of 
the  tip  vortex  for  this  test  condition  and  azimuth  as  shown  in  Figure 
24  to  be  0.015R  above  the  z/R  = -0.1  elevation.  If  it  were  shifted  by 
just  2.5  percent  of  the  radius  to  0.01R  below  the  measurement  line, 
good  correlation  with  the  radially  inward  velocity  would  result.  It  is 
thus  concluded  that  the  capability  of  the  Rotorcraft  Wake  Analysis  to 
predict  local  instantaneous  velocity  peaks  at  the  wake  boundary  depends 
upon  a very  accurate  positioning  of  the  tip  vortices  in  the  wake  repre- 
sentation. 


Regarding  the  lack  of  correlation  in  Figure  29  for  test  condition  1, 
the  test  data  indicating  a strong  peak  downwash  at  the  wake  boundary  may 
be  questioned  due  to  the  lack  of  a similar  peak  for  the  radial  component 
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for  this  condition,  and  thereby  questioning  the  presence  of  a close 
tip  vortex. 


In  summary,  analysis  of  the  comparisons  of  the  predicted  instan- 
taneous induced  velocities  with  the  Boatwright  test  data  revealed  the 
following: 

1.  The  correlation  of  the  generalized  wake  results  was  generally 
good. 

2.  Where  correlation  was  lacking,  it  could  generally  be  attributed 
to  known  shortcomings  of  the  generalized  wake  representation  (e.g., 
the  far  wake)  or  to  the  sensitivity  to  wake  excursions  during 
testing  due  to  ambient  wind  effects  and  the  0.5  deg  difference  in 
twist  between  the  two  blades.  However,  due  to  the  large  measured 
wake  excursions  and  the  uncertain  test  coning  angle,  the  direct 
correlation  between  the  wake  representation  generalized  from  model 
tests  and  the  wake  geometry  of  the  full-scale  0H-13E  rotor  could 
not  be  concluded. 

3.  The  use  of  the  generalized  wake  resulted  in  a considerable  im- 
provement in  the  correlation  over  that  of  the  undistorted  wake. 

Comparison  With  Test  Data  of  Hey son,  NACA  TN  D-393 


Induced  velocity  data  for  a two-bladed,  10-ft-diameter  rotor  model 
tested  on  a whirl  stand  mounted  in  the  full-scale  tunnel  at  the  NASA- 
Langley  Research  Center  are  presented  in  the  following  reference  (Ref- 
erence 19) ; 

Heyson,  H.H. , Measurements  of  the  Time-Averaged  and  Irstanta- 
neous  Induced  Velocities  in  the  Wake  of  a Helicopter  Rotor 
Hovering  at  High  Tip  Speeds,  NACA  Technical  Note  D-393 » July 
I960. 

This  report  presents  radial  distributions  of  the  time-averaged  component 
of  induced  velocity  normal  to  the  rotor  disk  at  two  elevations  below  the 
rotor  (z/R  = -0.145  and  -0.45).  Measurements  were  taken  in  and  out  of 
the  wake  of  the  hovering  rotor  using  rakes  of  total-head  and  static  pres- 
sure tubes.  In  addition,  a constant-temperature  hot-wire  anemometer  was 
installed  at  one  position  (z/R  = -0.1^5,  r/R  = 0.67)  to  record  the  time 
history  of  the  instantaneous  velocities  essentially  normal  to  the  rotor 
disk. 
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Information  on  the  rotor  design  and  test  conditions  are  presented 
in  Table  3- 


TABLE  3.  ROTOR  AND  TEST  INFORMATION  — HEYSON,  NACA  TN  D-393 


Rotor  Type 

Single,  model,  teetering  rotor 
Rotor  Design 


Number  of  blades 

2 

Rotor  solidity 

0.10 

Rotor  radius,  ft 

5.0 

Blade  twist,  deg 

-12.0 

Tip  chord,  ft 

0.63 

Airfoil  section 

NACA  0012 

Blade  taper  ratio 

2.0 

Selected  Test  Conditions 

Condition  1:  fi  = 0, 

fiR  = 700  fps. 

0?5  = ^ deg,  CT  = 

0.00225*,  v 

o 

= 23.5  fps 

Condition  2;  fi  = 0, 

flR  = 700  fps, 

#75  = 8 dee>  cT  = 

o.oo64*,  v = 
o 

39-5  fps 

Condition  3:  M = 0, 

flR  = 700  fps. 

e?5  = 12  deg,  CT  i 

= 0.0106*,  v 

o 

= 51.0  fps 

Condition  4:  m = 0, 

flR  = 500  fps, 

#75  = 8 dee>  CT  s 

0.0058*,  v = 
O 

27-0  fps 

* C not  indicated  in  reference:  calculated  values  listed. 
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Data  for  four  test  conditions  were  selected  which  include  variations  of 
collective  pitch  and  tip  speed.  Measured  values  of  thrust  coefficient 
were  not  included  in  the  »'eference.  The  calculated  values  using  the  test 
collective  pitch  with  both  the  (Generalized  and  classical  wake  representa- 
tions were  similar  (within  approximately  3 percent),  and  are  listed  in 
Table  3 alonf'  with  the  corresponding  momentum- induced  velocity  values. 

Comparisons  of  the  measured  and  theoretical  time-averaged  normal 
velocities  at  the  two  vertical  distances  (z/R  = -0.145  and  -0.45)  are 
presented  in  nondimensional  form  in  Figures  30  through  33.  In  Figures 
30  and  31j  comparisons  of  the  classical  and  generalized  wake  results 
with  the  test  data  are  presented  for  three  values  of  collective  pitch 
and  a constant  tip  speed  of  700  fps.  Although  the  magnitude  of  the 
theoretical  downwash  distribution  within  the  wake  is  generally  less 
than  measured,  the  general  characteristics  of  the  velocity  distributions 
and  the  influence  of  collective  pitch  (thrust  level)  are  predicted 
accurately.  As  expected,  relative  to  the  classical  wake  results  in 
Figure  30,  the  results  in  Figure  31  based  on  the  generalized  wake  are 
in  much  better  agreement  with  the  test  data  near  the  wake  boundary. 


The  disagreement  in  magnitude  of  the  downwash  based  on  the  general- 
ized wake  is  at  least  partially  attributable  to  the  lack  of  experimental 
wake  geometry  data  for  tapered  blades . Limited  wake  geometry  data  ac- 
quired in  the  test  reported  in  Reference  11  for  a model  rotor  with  six 
tapered  blades  (2:1  taper  ratio)  indicate  an  effect  of  taper  on  wake 
geometry,  but  the  data  were  not  of  sufficient  extent  to  be  generalized. 

The  capability  of  the  Rotorcraft  Wake  Analysis  to  predict  the 
insensitivity  of  nondimens ionalized  induced  velocity  to  tip  speed 
variations,  evident  in  the  test  results,  is  shown  in  Figures  32  and  33- 
In  these  figures,  results  for  test  conditions  2 and  4 are  compared, 
which  correspond  to  tip  speed  variations  between  700  and  500  fps, 
respectively,  a collective  pitch  value  of  8 deg,  and  similar  thrust 
coefficients.  Independence  of  tip  speed  at  similar  thrust  coefficient 
levels  is  predicted  by  both  the  classical  and  generalized  wake  results. 
These  results  are  consistent  with  the  statement  of  Reference  19  and  the 
experimental  wake  generalization  of  Reference  8 that  tip  speed  ( at 
least  up  to  700  fps)  has  essentially  no  independent  effect  on  the  wake 
of  the  hovering  rotor. 

In  Figure  34,  the  time  history  of  the  instantaneous  induced  velocity 
component  normal  to  the  rotor  disk  measured  at  one  point  in  the  rotor 
wake  with  a hot  wire  is  reproduced  directly  from  the  reference.  Although 
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no  scale  is  included,  the  shape  of  the  2/rev  velocity  pulses  associated 
with  the  two-bladed  rotor  is  clearly  defined.  For  comparison  of  the 
shape  of  the  velocity  pulses,  the  velocity  time  history  predicted  using 
the  generalized  wake  was  normalized  and  phased  to  the  test  data.  As 
shown  in  Figure  34,  the  shapes  of  the  theoretical  and  experimental  time 
histories  are  similar.  Also  presented  in  Reference  19  are  the  recorded 
root-mean-square  values  of  the  velocity  pulses.  For  the  test  condition 
shown  (Condition  3),  the  velocity  fluctuations  are  recorded  as  13  per- 
cent of  the  local  mean  induced  velocity  and  25  percent  of  the  momentum 
induced  velocity.  The  predicted  fluctuations  were  found  to  be  in  favor- 
able agreement  with  these  measured  values,  in  that  the  respective  pre- 
dicted fluctuations 


1/2  v 


1/2  V -V, 
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were  calculated  to  be  15  percent  and  26  percent. 

Comparison  With  Test  Data  of  McKee  and  Naeseth,  NACA  TN  4239 

Induced  velocity  data  in  the  form  of  local  dynamic  pressure  measure- 
ments for  a two-bladed,  6-foot-diameter  rotor  model  with  constant  chord, 
untwisted  blades  are  presented  in  the  following  reference  (Reference  20); 

McKee,  J.W.  and  R.L.  Naeseth;  Experimental  Investigation  of  the 
Drag  of  Flat  Plates  and  Cylinders  in  the  Slipstream  of  a Hover- 
ing Rotor,  NACA  TN  4239  > April  1958. 

The  model  rotor  was  tested  in  a partially  enclosed  outdoor  area.  Dynamic 
pressure  data  were  taken  using  both  Pitot  tubes  and  small  disks  at  var- 
ious radial  and  vertical  locations  on  the  slipstream  side  of  the  hovering 
rotor.  Local  dynamic  pressures  normal  to  the  rotor  disk  were  measured 
which  are  proportional  to  the  time-averaged  normal  component  of  induced 
velocity  squared. 

Information  on  the  rotor  design  and  single  test  condition  is  pre- 
sented in  Table  4.  The  test  was  run  at  a tip  speed  of  367  fps  and  a 
collective  pitch  of  11  deg,  and  rotor  performance  was  measured. 

Sample  dynamic  pressure  test  data,  reproduced  directly  from  the 
reference,  are  presented  in  Figure  35-  Differences  in  the  data  between 
the  two  measurement  devices  and  data  fluctuations  excited  by 
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recirculation  from  the  surrounding  walls  and  ambient  wind  are  cited 
in  the  reference.  Random  variations  of  the  time  histories  were  espe- 
cially large  near  the  wake  boundary,  and  it  was  consistently  noted 
that  the  slipstream  was  not  centered  with  respect  to  the  rotor. 


TABLE  4.  ROTOR  AND 

TEST  INFORMATION  — McKEE  & NAESETH, 

NACA  TN  4239 

Rotor  Type 

Single,  model,  teetering  rotor 

Rotor  Design 

Number  of  blades 

2 

Blade  taper  ratio 

1.0 

Rotor  radius,  ft 

3-0 

Blade  aspect  ratio 

9-1 

Blade  chord,  ft 

0.33 

Blade  twist,  deg 

0.0 

Rotor  solidity 

0.071 

Airfoil  section 

NACA  0012 

1 

I Test  Condition 

M = 0,  flR  = 367  fps. 

fl75  * U 

.0  deg,  C = 0.00602,  v = 
I O 

20.13  fps 

4 


Comparisons  of  the  predicted  and  measured  dynamic  pressures  are 
presented  in  Figures  35  and  36  for  five  elevations  ranging  from  z/R  = 
-0.104  to  -0.993.  Figure  36  includes  the  data  variations  from  the  two 
measurement  devices  and  the  unsymmetrical  side-to-side  variations  due 
to  wake  unsteadiness.  The  predicted  and  test  results  are  generally  in 
good  agreement  near  the  rotor.  The  magnitudes  of  the  predicted  results 
are  generally  somewhat  lower  than  the  test  results,  although  the 
dynamic  pressure  differences  relative  to  differences  in  induced  veloc- 
ities are  accentuated.  For  example,  a 20  percent  difference  in  dynamic 
pressure  is  representative  of  a 10  percent  difference  in  induced  velo- 
city due  to  the  square  relation.  Also,  to  achieve  the  measured  thrust 
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coefficient  in  the  Rotorcraft  Wake  Analysis,  the  collective  pitch  was 
lowered  from  the  11  deg  stated  in  the  reference  to  10.4  deg.  Ii  ‘'he 
11  deg  had  been  used  in  the  analysis,  the  predicted  dynamic  pressure 
level  would  have  been  increased  and  would  have  agreed  more  closely 
with  the  test  data.  However,  assuming  the  test  information  to  be 
accurate,  the  analytical  inconsistency  is  recognized,  and  a preferred 
choice  between  the  two  approaches  is  not  apparent.  The  generalized 
wake  results,  which  properly  account  for  near-wake  contraction,  are 
in  better  agreement  with  the  classical  wake  results  near  the  rotor. 

The  unsteady  behavior  and  vortex  core  dissipation  in  the  real  hovering 
far  wake,  not  represented  in  the  analytical  wake  models,  are  evident  in 
the  test  data  by  the  increased  region  of  test  data  and  lower  and 
broadened  pressure  peaks  at  z/R  = -0.660  and  -0.993* 


Comparison  With  Test  Data  of  Fradenburgh 


Experimental  induced  velocity  contours  for  a two-bladed,  2-foot 
diameter  rotor  model  tested  on  an  indoor  hover  test  stand  are  presented 
in  the  following  reference  (Reference  21); 

Fradenburgh,  E.  A.,  Flew  Field  Measurements  for  a Hovering 
Rotor  Near  the  Ground.  Proceedings  of  the  5th  Annual  Western 
Forum  of  the  American  Helicopter  Society,  Los  Angeles,  California, 
September  1958. 

Although  the  data  from  this  reference  have  some  limitations  for  correla- 
tion purposes,  they  were  selected  recognizing  that  the  model  test  was  one 
of  the  first  to  provide  a comprehensive  set  of  downwash  data  beneath  a 
hovering  rotor  and  has  remained  over  the  years  as  one  of  the  few  pub- 
lished sources  of  such  data.  The  limitations  are  associated  with  the 
small  rotor  size,  low  aspect  ratio,  and  the  use  of  Pitot  tubes  to  mea- 
sure the  time-averaged  velocity  component  normal  to  the  rotor  disk  in  a 
fluctuating  flow  field. 

Rotor  and  test  information  are  presented  in  Table  5*  The  rotor  was 
tested  out-of-ground-effect  at  one  test  condition  for  which  the  thrust 
coefficient  was  approximately  0.007*  Collective  pitch  (675)  was  not  in- 
dicated in  the  reference.  The  9 75  values  used  in  the  analysis  were  the 
values  required  to  produce  the  test  thrust  coefficients  (0*^  =»  9*2  deg 
for  classical  wake;  Q7C-  =»  10.0  deg  for  generalized  wake). 


The  experimental  downwash  velocity  contours  measured  at  four  loca- 
tions below  the  rotor  (indicated  by  the  dotted  lines),  and  nondimension- 
alized  by  the  momentum  induced  velocity  (v  =»  35*5  fps),  are  presented  in 
Figure  37* 
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TABLE  5. 


ROTOR  AND  TEST  INFORMATION  — FRADENBURGH 


1 


Rotor  Type 

Small,  single,  model, 

hingeless 

rotor 

Rotor  Design 

Number  of  blades 

2 

Blade  taper  ratio 

1.0 

Rotor  radius,  ft 

1.0 

Blade  aspect  ratio 

6.0 

Blade  chord,  ft 

0.167 

Blade  twist,  deg 

0.0 

Rotor  solidity 

0.106 

Airfoil  section 

NACA  0012 

Test  Condition 

ix  = 0,  ftR  = 600  fps 

> #75  not 

indicated,  C^  s 0.007, 

v =35-5  fps 

O 

Velocity  contours  calculated  using  the  classical  and  generalized  wake 
representations  in  the  Rotorcraft  Wake  Analysis  are  presented  in 
Figures  38  and  39-  Direct  comparisons  of  the  analytical  time-averaged 
velocities  with  the  test  velocities  can  be  made  by  comparing  the  con- 
tours shown  in  Figures  38  and  39  with  those  of  Figure  37. 


Comparisons  of  the  predicted  velocities  with  the  test  data  reveal 
the  following: 

1.  The  general  character  of  the  analytical  and  experimental  veloc*- 
ity  contour  lines  are  similar. 


| 

[ 


2.  As  anticipated,  the  velocity  contours  based  on  the  classical 
wake  geometry  are  shifted  outward  from  the  experimental  contours 
below  the  blade  tip  region  due  to  the  absence  of  wake  contraction. 
It  is  noted  that  the  experimental  velocities  do  not  show  evidence 
of  the  expansion  of  the  boundaries  of  the  far -wake  region  evident 
in  other  sources  of  data.  (A  wake  photograph  in  Reference  21 
indicates  expansion  at  z/R  = -1.2.) 


I 
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3.  On  e.  point-to-point  basis,  the  analytical  velocity  contours 
based  on  the  generalized  wake  geometry  are  consistently  lower  than 
the  experimental  values  by  an  increment  of  approximately  0.5  vQ. 


I 


4.  A larger  region  of  upflow  beneath  the  rotor  hub  is  noted  for 
the  analytical  results  which  is  attributable  to  the  establishment 
of  the  inboard  wake  boundary  in  the  analysis  below  the  root  cutout 
location,  and  to  the  lack  of  a provision  for  mixing  of  the  flow  in 
this  low-pressure  region. 

There  are  several  reasons  for  the  third  observation,  and  they  are 
primarily  related  to  the  accuracy  of  the  analytical  wake  representation. 

A comparison  of  the  generalized  wake  representation  used  for  this  rotor 
in  the  analysis  with  the  wake  geometry  indicated  in  a photograph  in  the 
Fradenburgh  reference  rp;eaiod  significant  differences.  This  is  believed 
to  be  attributable  to  the  extrapolation  of  the  generalized  wake  results 
of  Reference  8 acquired  for  two-bladed  rotors  with  relatively  low  solidi- 
ties (0.035  and  0.047),  resulting  from  high  aspect  ratios  (18.2  and  13.6), 
to  a rotor  with  a solidity  of  0.1  and  u.n  aspect  ratio  of  6.  Recent  indi- 
cations are  that  the  generalized  wake  equations  (Appendix  C)  do  not  apply 
to  rotors  with  low-aspect-ratio,  high-solidity  blades.  For  example, 
indications  a^e  that  the  tip  vortex  axial  coordinates  (k^  in  Appendix  C) 
should  be  generalized  on  a parameter  such  as  CT/b  rather  than  CT/a,  to 
apply  to  such  rotors. 

Comparison  With  Test  Data  of  Flemming 

The  previously  presented  sources  of  hovering  rotor  data  have  all 
been  for  rotors  with  two  blades.  In  order  to  compare  the  induced  veloc- 
ities from  th2  Rotorcraft  Wake  Analysis  with  test  data  for  a large  number 
of  blades,  hot-film  anemometer  data  for  a seven-bladed  rotor  were  selec- 
ted from  Reference  22. 

Flemming,  R.  J.,  l/20  Scale  Model  Rotor  Wake  Survey,  Internal 
Sikorsky  Documents,  May  1969  and  March  1972. 


The  Sikorsky- sponsored  test  program  was  conducted  at  the  UTRC  model  rotor 
hover  test  facility.  Hot  film  anemometer  equipment  was  employed  to  mea- 
sure the  time-averaged  downwash  velocity  at  several  collective  pitch 
angles  (thrust  levels)  and  two  tip  speeds.  Traverses  of  a hot-film  probe 
were  made  at  several  elevations  below  the  rotor  to  result  in  radial  dis- 
tributions of  the  induced  velocity  component  normal  to  the  rotor  disk. 
Information  on  the  rotor  design  and  four  selected  test  conditions  is  pre- 
sented in  Table  6. 


u --A:  L/- 


TABLE  6 . - ROTOR  AND  TEST  INFORMATION  — FLOWING 


Rotor  Type 

Single,  model,  semi articulated  (flapping)  rotor 


Potor  Design 
Number  of  blades 
Rotor  radius,  ft 
Blade  chord,  ft 
Rotor  Solidity 


7 Blade  taper  ratio 

1.975  Blade  aspect  ratio 

0.108  Blade  twist,  deg 

0.122  Airfoil  section 


NACA  0012 


Test  Conditions 

Condition  1:  ft  = 0,  HR  = 620  fps,  = 9 deg,  CT  = O.OO69,  vq  = 36.3  fps 

Condition  2:  ft  = 0,  OR  = 620  fP3J  $75  = 10  de8>  CT  = 0.0079,  Vq  = 39-0  fps 

Condition  3:  ft  = 0,  flR  - 620  fps,  6^  = 12  deg,  CT  = 0.0103,  vq  = 44.5  fps 

Condition  4:  ft  =»  0,  flR  = 465  fps,  6^  = 9 deg,  CT  = O.OO67,  vq  = 27-0  fps 


The  Rotorcraft  Wake  Analysis  was  run  with  classical  and  generalized 
wake  models  for  the  four  test  conditions.  When  necessary,  collective 
pitch  was  adjusted  in  the  analysis  from  the  test  values  to  produce  the 
measured  thrust  coefficients.  Comparisons  of  the  theoretical  and 
measured  time-averaged  induced  velocities  for  each  test  condition  at 
two  distances  below  the  rotor  hub  center  (z/R  = -0.158  and  -O.316)  are 
presented  in  Figures  40  to  43.  Similar  to  previous  comparisons  with 
other  data  sources,  both  wake  models  in  the  analysis  correctly  predict 
the  character  of  the  radial  distributions  of  velocity.  The  analytical 
results  in  the  outboard  region  reflect  the  lack  of  contraction  of  the 
classical  wake  and  the  smooth  contraction  to  O.78R  of  the  generalized 
wake.  The  general  increase  in  inboard  flow  velocity,  at  fixed  radial 
coordinates,  with  increased  distance  from  the  rotor  is  consistently 
predicted  for  this  rotor.  Figure  44  is  included  to  show  that  this  applies 
over  the  full  range  of  z/r  elevations  (z/r  = -O.O85  to  -0.99)  for  which 
measurements  were  made.  For  clarity,  results  for  three  intermediate  posi- 
tions (z/R  = -0.158  and  -0.3l6  presented  in  Figure  4l  and  -0.470)  have 
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been  omitt'-d,  but  t.h«-  analyt ical  v* •locitif:  in  the  lip.-troam  progren- 
. i v ly  inert •?..•••<!  in  magni  tud*-  with  increaning,  di  "lance  from  the  rotor 
con.i.ten4  with  the  t>  nt  renultn. 


of  the  generalized  wake  representation  ov< -rprodict the  absolute 
mafnituie  of  th>-  induced  velocity,  and  the  clarrical  wake  result."  gener- 
ally indi cal  ••  better  correlation  inboard  than  the  generalized  wake 
re.ul+.  for  thin  rotor.  Howevc-r,  due  to  the  large  number  of  waken  for 
.■eV'-n  blade.-,  and  the-  clone  pannage  of  the  tip  vortices  to  the  blade  , the 
r*-. -\.lt  r wer--  found  to  be  very  sensitive  to  wake  repre  rent  at  i on  (tip  vortex 
geometry,  analytical  modeling  of  tip  vort  --x  roll-up,  and  selection  of 
radial  point..-  of  origination  of  tip  vortex  filament,  from  the  blade.-). 
Also,  t.  h * time -averaging  of  indue-  1 v-  locit  i--n  calculat'  d at  discrete  time 
incr--m--n4  • in  -n.it  iv  'o  tie-  -l--ction  of  blad'-  azimuth  increment  for 
+ he  clo  ■ ly  . paced  .-•-•v< -n  fcla  i-  . Tn  addition,  th<-  clone  proximity  of  nuc- 
c-  . .'iv-  tip  vortices  for  a . -v-  n-bl  ad • ■ 1 rotor  r-ultn  in  tip  vortex  mixing, 
wak  • ur.  *- • -adi n-  • , and  4 ip  vor*  • x cor--  di  ipat.ion  cloner  to  the  rotor 
4 bar.  for  rdor.-  with  lew-  r numb-  rn  of  bind--.  . Thi . trend  with  incr-  a.  ing 
m.mb-r  of  blab-'  ; . -vid-nt  in  t.h"  flow  visualization  r*  nultn  shewn  in 
f-  r-  t.c-  in  which  th-  wak--.  of  two-,  four-,  nix-  and  ■ ipht-bladed 
ro+or  nr--  pi--  -'-nt-  d,  "or  -■xnmpl-.-,  t).-  tip  vorti  c-  n of  ti:e  eight-bladed 
.'••o'- or  v-  found  *o  fc- com-  ur.  4 ■ ady  at  about  0.2R  from  the  rotor  compared 

4o  C .ij  4o  '.5°  for  th-  4 wc-blad-d  rotor,  in  gn-ral,  th-.-  results  indicate 
4 h-  importanc--  of  accurat-ly  moi  ling  th'-  complex  wak  - of  hovering  rotor.: 
with  larr--  r.umb--r  • of  bind-  . . 


Imparl  son  '-'it:  ;'-.(.-r  .(.-c  i a1  a < f Miller,  .ass.  a;..1  Perlmutter. 


n <-r -or  •.’.  compare  tne'.  re-  i cal  isiuee:  olocivios  wl’.i;  data  from 

’so  fi  i ’ os'  ' f a full -scale-  ueli  <:• pter,  ■ to  following  reference 

fei'e.  Ce  sr.]_,  C*  0 U 


Mi  Her , a:.'  an  i r'-rlir.u' *.er.  ."heorotical  ar.d 

experimental  :.  . es-.i  -at ion  of  tig.-  Ins*  antaneous  induced  Veloci 
:-ioli  i:.  the  V.'ake  of  a Lifting  'rotor,  gyr.asciences  L’orp..  VCAAVLA'-v 
ee.-jsiual  p. -pert  • 7-'- , i-ustis  directorate'.  Army  Air  debility 

•enearen  a.-.:  : eyeiopiricr.t  Las oratory . lor'  ustis.  7a. . January  19'' 
A!  • 


Veloci'...'  ’ect  data  were  measured  in  flig-g  on  an  Army  i 7-55A  helicopter 
using  : cti.  pressure  (Pitot  pro:  os)  ar.  i ho*  - film  instrumentation  mounted 
on  a radially  oxter,  iir.g  : ton  eo-a‘.::  the  ro'.or  at  various  azimuth 
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positions  and  axial  distances  from  the  rotor  disk.  Both  time-averaged 
and  instantaneous  velocity  data  were  acquired.  In  the  course  of  con- 
ducting this  correlation  study,  it  was  found  that  the  data  in  this  ref- 
erence was  not  very  useful  due  to  the  imprecise  manner  of  data  presenta- 
tion, the  limited  extent  of  the  data,  and  the  degree  of  scatter  in  some 
of  the  data.  In  particular,  it  was  found  that  the  instantaneous  veloc- 
ity data  were  presented  in  unusable  form. 

Rotor  and  hover  flight  test  information  is  presented  in  Table  7* 
Collective  pitch  was  not  documented  in  the  reference,  and  was  thus  es- 
tablished in  the  analysis  as  6.8  deg  to  produce  a thrust  coefficient 
consistent  with  the  nominal  gross  weight  of  the  aircraft. 


TABLE  7-  ROTOR  AND  TEST  INFORMATION  — 
MILLER,  TANG  & PERIMUTTER,  USAAVLABS  TR  67-68 


Rotor  Type 

TH-55A,  full-scale,  articulated,  rotor 
Rotor  Design 

Number  of  blades  3 

Rotor  radius,  ft  12.6U5 

Blade  chord,  ft  0.80 

Rotor  solidity  0.060 

Hover  Test  Condition  (Flight  Test) 

/a  = 0,  UR  = 596  fps,  9„r-  not  indicated,  C = 0.00377*,  v = 25.9  fps 

75  10 


Blade  taper  ratio 
Blade  aspect  ratio 
Blade  twist,  deg 
Airfoil  section 


1.0 

15.9 

-8.0 

NACA  0015 


* Based  on  nominal  gross  weight. 


The  three  components  of  predicted  time-averaged  induced  velocity  are 
compared  with  both  the  hot-film  and  pressure  data  for  the  hover  condition 
in  Figures  45  and  46.  The  velocities,  nondimensionalized  by  the  momentum 
value,  are  presented  for  the  measurement  radial  stations  between  0.55R 
and  0.80R  at  two  distances  below  the  rotor  (z/R  = -0.26  and  -0.45).  The 
theoretical  results  are  shown  to  correlate  well  within  the  scatter  and 
limited  extent  of  the  test  data.  The  limited  radial  extent  of  the  test 
data  precludes  the  opportunity  to  comparatively  assess  the  results  based 
on  the  classical  and  generalized  vake  representations  at  and  outboard  of 
the  wake  boundary. 


CORRELATION  STUDY  FOR  FORWARD- FLIGHT  CONDITIONS 


Five  of  the  selected  sources  of  data  (References  23  to  27)  listed 
in  Table  1 contain  information  pertaining  to  forward-flight  test  condi- 
tions. Comparisons  of  results  of  the  Rotorcraft  Wake  Analysis  with 
test  data  from  each  are  presented  below. 


Comparison  With  Test  Data  of  Hey son  and  Katzoff,  NACA  TR  131S 


Induced  velocity  data  for  a 15-foot  diameter,  two-bladed  model 
rotor  tested  in  the  NASA  Langley  full-scale  wind  tunnel  are  presented  in 
the  following  reference  (Reference  24): 


Heyson,  H.  H.  and  3.  Katzoff,  Induced  Velocities  Near  a Lifting  Rotor 
With  Nonuniform  Disk  Loading,  NACA  Technical  Report  1319,  1957- 


Presented  in  this  reference  are  time-averaged  values  of  the  component  of 
induced  velocity  normal  to  the  rotor  disk.  These  measurements  were 
acquired  by  surveying  a large  number  of  points  above,  below,  and  aft  of 
the  rotor  with  a rake  of  pressure  tubes. 


Information  on  the  rotor  design  and  test  conditions  is  presented 
in  Table  8,  As  noted  in  the  table,  data  for  three  test  conditions  with 
advance  ratios  of  0.095,  0.l4o  and  0.232  were  selected  for  this  correla- 
tion study.  An  undistorted  waRe  model  with  a rolled-up  tip  vortex  was 
used  in  the  calculations.  Sample  results  of  the  correlation  study  are 
presented  in  Figures  47  through  53.  These  figures  contain  comparisons 
of  measured  and  theoretical  values  of  the  time-averaged  induced  velocity 
ratio  (normal  to  the  rotor  disk)  at  selected  field  points.* 


*The  comparison  is  made  with  test  data  exactly  as  presented  in  the  refer- 
ence, in  which  downflow  is  positive.  To  be  consistent  with  the  sign 
convention  used  herein  for  vz/v0  (downflow  negative),  the  notation 
-vz/v0  is  used  in  the  figures. 
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TABLE  8. 


ROTOR  AND  TEST  INFORMATION  — HEYSON  & KATZOFF,  NACA  TR  1319 


Rotor  Type 

Single,  model,  teetering  rotor 


Rotor  Design 
Number  of  blades 
Rotor  radius , ft 
Blade  chord,  ft 
Rotor  solidity 

Test  Conditions 


2 

7-5 

0.64 

0.0543 


Blade  taper  ratio 
Blade  aspect  ratio 
Blade  Twist,  deg 
Airfoil  section 


1 

11.7 

0 

NACA  0012 


Condition  1:  /z  = 0.095*  ftR  = 5°0  fps,  C^,  = 0.00320,  a = 

Condition  2:  /*  = 0.140,  ftR  = 500  fps,  C^  = 0.00371*  « = 

Condition  3 : /x  = 0.232,  ftR  = 450  fps,  CT  = 0.00321,  a = 

Blade  control  angles  not  indicated  in  reference. 


-9.?°,  vc 
-5.3°,  vc 
-9.5°,  v 


8.0  fps 

6.50  fps 
3.04  fps 


For  Figures  47  through  49  the  advance  ratio  is  0.140,  the  thrust 
coefficient  is  0.0037,  and  the  wake  skew  angle  is  82.3  degrees  based  on 
a momentum  induced  velocity  of  6.5  fps.  In  Figure  47,  the  induced  veloc- 
ity ratios  are  compared  at  several  distances  above  and  below  the  rotor 
disk  in  the  lateral  plane  of  symmetry  of  the  rotor.  In  Figure  48,  com- 
parisons are  made  at  several  distances  fore  and  aft  of  the  rotor  center 
and  in  a plane  parallel  to  and  0.07R  beneath  the  rotor.  This  plane  was 
selected  for  comparison  on  the  basis  that  it  contains  the  measured  points 
that  are  closest  to  the  rotor  and  within  the  rotor  wake.  The  trends  of 
the  velocity  distributions  are  predicted  accurately,  and  good  correlation 
is  noted  except  at  some  points  beneath  the  advancing  side  of  the  rotor. 

In  Figure  49,  the  predicted  values  are  compared  with  the  data  as  inter- 
polated in  NACA  TR  1319  along  the  longitudinal  and  lateral  centerlines  of 
the  rotor  disk.  Good  correlation  is  noted  along  the  longitudinal  center- 
line  except  at  x/R  = 3 *14,  for  which  the  number  of  wake  revolutions  used 
in  the  Rotorcrafs  Wake  Analysis  was  insufficient  (four  wake  revs  were  used 
which  truncated  the  wake  in  the  vicinity  of  x/R  = 3»l4).  However,  lack  of 
correlation  in  this  downstream  region  is  also  expected  due  to  the  down- 
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stream  distortions  of  the  wake.  Reasonable  correlation  is  noted  on  the 
retreating  side  of  the  lateral  centerline.  However,  on  the  advancing 
side,  the  induced  velocity  is  overpredicted  as  is  consistent  with  the 
previous  figures. 

In  Figure  50,  the  induced  velocity  ratio  is  compared  for  the  test 
condition  defined  by  an  advance  ratio  of  0.095  and  a thrust  coefficient 
of  0.0032.  The  undistorted  wake  skew  angle  for  this  condition  is  75 
deg,  based  on  a momentum  induced  velocity  of  8.0  fps.  The  degree  of  cor- 
relation is  generally  similar  to  that  at  the  0,l40  advance  ratio. 
Comparisons  of  measured  and  predicted  velocities  for  the  0.232  advance 
ratio  condition  are  presented  in  Figures  51  to  53.  The  undistorted 
wake  skew  angle  for  this  condition  is  84  degrees  and  the  momentum  induced 
velocity  is  3.0  fps.  Since  at  high  rotor  advance  ratios  the  wake  is 
swept  back  rapidly,  leaving  the  forward  half  of  the  rotor  in  relatively 
clean  flow  (few  blade-vortex  interactions),  comparisons  were  made  at 
points  in  the  aft  region  of  the  rotor.  Comparisons  are  presented  at  the 
lateral  plane  of  symmetry,  x/R  = 0,  and  at  the  x/R  = 0.5  and  1.07  planes. 
It  is  noted  that  the  experimental  accuracy  for  this  condition  was  recorded 
in  the  reference  as  a -25  percent  compared  to  the  110  percent  for  the 
other  conditions. 

Generally,  the  degree  of  correlation  appears  to  be  reasonable,  con- 
sidering the  undistorted  wake  representation  used  in  the  analysis. 
Significant  distortions  and  mixing  of  the  wake  on  the  lateral  portions 
of  the  rotor  (roll-up  into  a pair  of  vortices  similar  to  a fixed  wing) 
are  evident  from  the  velocity  contours  presented  in  the  reference  and 
more  recent  flow  visualization  results.  Whether  the  discrepancies  that 
do  exist  (particularly  on  the  advancing  side  at  x/R  = 0,  z/R  = 0 and 
-0.07)  s.re  due  to  wake  distortions,  the  finite  filament  wake  represen- 
tation, inaccuracies  in  the  input  analytically  determined  controls  and 
flapping  that  affect  lateral  trim,  or  sensitivity  to  other  assumptions 
or  limitations  of  the  analysis,  remains  to  be  determined. 

Comparison  With  Forward  Flight  Test  Data  of  Miller,  Tang,  and  Perlmutter, 
USAAVLABS  TR  67-68 

In  addition  to  the  comparison  with  the  hover  test  data,  presented 
previously,  induced  velocities  from  the  Rotorcraft  Wake  Analysis  were 
compared  with  the  forward  flight  test  data  from  the  following  reference 
(Reference  23). 

Miller,  11.,  J.  C.  Tang  and  A.  A.  Perlmutter,  Theoretical  and 
Experimental  Investigation  of  the  Instantaneous  Induced  Velocity 
Field  in  the  Wake  of  a Lifting  Rotor,  USAAVLABS  Technical  Report 
67-^8,  January  19^8. 
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Rotor  and  forward  flight  test  information  are  presented  in  Table  9. 

As  noted  in  this  table,  the  two  flight  conditions  have  low  advance 
ratios,  O.OU9  and  0.11,  corresponding  to  flight  speeds  of  approximately 
17  and  4o  knots,  respectively. 


Radial  distributions  of  induced  velocity  were  predicted  for  the 
following  test  azimuth  and  axial  positions: 

Hot-Film  Boom  Azimuth  Positions,  deg  30,  90,  150,  210,  270,  330 
Axial  Position  (from  hub),  z/r  -0.22,  -0.26,  -0.49 

The  momentum  skew  angles  of  the  classical  wake  representation  were  X = 58 
and  80  deg,  corresponding  to  the  momentum  induced  velocity  values  of  18.6 
and  9.3  fps  for  the  17-  and  40-kt  conditions.  Results  for  the  two  flight 


TABLE  9.  ROTOR  AND  TEST  INFORMATION  — MILLER,  TANG  & PERLMUTTER, 
USAAVLABS  TR  67-68,  FORWARD  FLIGHT 


TH-55A,  full-scale,  articulated  rotor 


Rotor  Design 
Number  of  blades 
Rotor  radius,  ft 
Blade  chord,  ft 
Rotor  solidity 


3 

Blade  taper  ratio 

1.0 

12.645 

Blade  aspect  ratio 

15.8 

0.80 

Blade  twist,  deg 

-8.0 

0.060 

Airfoil  section 

NACA 

Forward  Flight  Test  Conditions  (Flight  Test) 

Condition  1:  fi  = 0.049,  QR  = 596  fps,  CT  = 0.00377*,  « = -0.2°,  vq  = 18.6  fps 
Condition  2:  fi  = 0.11,  flR  = 596  fps,  CT  = 0.00377*,  a = -2.3°,  vo  = 9-3  fpa 

Blade  control  angles  not  indicated  in  reference. 

* C^  based  on  nominal  gross  weight. 


-r-'i  •* & "/.Ckr*? ** rtJ* * Vi  f .>* ' *»  * i • « 


condition:-  are  compared  with  the  time-averaged  test  data  (normal  compo- 
nent) in  Figures  54  and  55.  (instantaneous  velocity  data  in  the  reference 
was  presented  in  unusable  form.)  As  revealed  in  these  figures,  the  degree 
of  correlation  varies  from  good  to  poor.  This  is  believed  to  be  partially 
attributable  to  the  increased  sensitivity  of  local  velocities  to  wake  dis- 
tortions for  low- speed  conditions.  To  illustrate  this,  two  plots  were 
selected  from  Figures  54  and  55  and  are  presented  in  Figure  56,  they  are 
representative  of  the  range  of  correlation.  The  test  results  at  the 
selected  90-degree  azimuth  ire  fairly  similar  at  the  two  flight  speeds. 
Although  the  theoretical  results  correlate  well  at  the  4o-kt  condition, 
a significant  difference  is  noted  at  the  two  near  positions  to  the  rotor 
for  the  17.4-kt  condition.  The  charge  from  large  downflow  values  to 
upflow  in  the  predicted  radial  distributions  for  these  two  positions  is 
indicative  of  the  close  proximity  to  and  crossing  of  the  wake  boundary. 

The  high  inboard  downwash  velocities  predicted  for  these  two  axial  posi- 
tions are  indicative  of  placement  within  the  wake.  The  velocity 
gradient  indicates  the  wake  boundary,  and  the  upflow  predicted  outboard 
indicates  placement  outside  the  wake.  For  the  furthest  axial  position 
(z/P  = -0.49),  the  velocity  distribution  changes  character  since  the 
points  beyond  approximately  0.6r  lie  outside  the  wake.  The  discrepancy 
is  thu.  believed  to  be  mainly  attributable  to  the  classical  wake 
representation  used  in  the  analysis,  for  flow  visualization  studies 
have  indicated  that  distortions  of  the  wake  boundary  are  most  severe  at 
low- speed  conditions.  Although  some  differences  are  expected  due  to 
uncertainties  in  the  rotor  test  condition  parameters  and  other  aceump- 
tions  in  the  analysis,  results  for  points  at  other  azimuth  positions 
appear  support  the  conclusion  that  the  degree  of  correlation  depends 
on  the  degree  of  wake  distortion.  In  general,  points  with  low  theoretical 
induced  velocity  magnitudes  are  outside  the  classical  wake,  and  points 
with  the  higher  magnitudes  are  within  the  wake. 

Comparison  With  Test  Data  of  Landgrebo  and  Johnson 

The  flow  velocity  data  of  primary  interest  to  the  rotor  aerody- 
namicist  are  the  instantaneous  velocities  in  the  vicinity  of  the  rotor 
blades.  The  accurate  measurement  of  these  velocities  has  until  recently 
presented  the  rotor  experimentalist  with  a dilemma  --  how  to  measure 
the  flow  at  the  rotor  disk  without  influencing  the  flow  and  without 
requiring  concern  for  the  physical  impact  of  the  blades  with  the  measure- 
ment device.  In  the  last  3 years,  it  has  been  demonstrated  that  the 
dilemma  can  be  solved  by  using  a laser  velocimeter  (LV).  The  experi- 
mental velocities  of  the  correlation  study  presented  below,  and  the 
one  following,  were  acquired  from  LV  measurements. 
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In  1973,  LV  techniques  were  applied  at  UTRC  to  a model  rotor  in  a 
wind  tunnel,  and  velocity  time-histories  were  measured  at  several  fixed 
points  in  close  proximity  to  passing  blades.  This  demonstration  test  is 
described  in  the  following  reference  (Reference  25). 

Landgrede,  A.  J.,  and  B.  V.  Johnson,  Measurement  of  Model  Heli- 
copter Rotor  Flow  Velocities  with  a Laser  Doppler  Velocimeter,  Jl. 
American  Helicopter  Society,  Vol.  19,  No.  3,  July  1974. 


A photograph  of  the  rotor  and  LV  system  in  the  UTRC  4-  x 6-foot  wind 
tunnel  is  shown  in  Figure  57-  The  velocity  measurement  point  is  the 
focal  point  (intersection)  of  the  two  laser  beams  shown.  Rotor  design 
and  test  information  are  presented  in  Table  10.  A schematic  showing 
the  side  view  of  the  two-bladed  rotor  with  the  blades  in  the  fore  and 
aft  positions  ($=  0 and  180  deg)  for  the  0.15  advance  ratio  condition  is 
presented  in  Figure  58.  The  test  was  conducted  as  a demonstration  to 
show  the  feasibility  of  applying  LV  techniques  in  a wind  tunnel  to  meas- 
ure rotor  instantaneous  flow  velocities,  and  thus  only  a limited  number 
of  velocity  measurements  were  taken.  Velocity  data  were  measured  at 


TABLE  10.  ROTOR 

AND  TEST 

INFORMATION  --  LANDGREBE 

& JOHNSON 

Rotor  Type 

Single,  model,  articulated  rotor 

Rotor  Design 

Number  of  blades 

2 

Blade  taper  ratio 

1 

Rotor  radius,  ft 

2.15 

Blade  aspect  ratio 

17.6 

Blade  chord,  ft 

0.122 

Blade  twist,  deg 

0 

Rotor  solidity 

0.0363 

Airfoil  section 

NACA  0012 

Test  Condition 

fi  = 0.15,  QR  = 300  fps 

, cT  = 0.0026,  0?5  = 9°,  = 2.5 

fps, 

0 0 

,0 

_ 0 „ n° 

Ai  = 0 > B!  - 0 . a0 

s s 

= 3 , ax 

s 

= 3.2  , n = -12.8 

points  in  the  rotor  vertical  plane  of  symmetry  forward  of  the  rotor  hub 
for  several  vertical  distances  from  the  rotor  disk.  These  vertical  dis- 
tances ranged  from  0.08R  above  the  disk  to  0.16R  below  the  disk  (mea- 
sured from  the  flapped  blade  at  $ = 180  deg),  as  shown  in  Figure  58. 
Vertical,  vz<  , and  streamwise,  vxt,  velocity  components  were  measured. 

To  assist  in  the  interpretation  of  data,  flow  visualization  photo- 
graphs (similar  to  Figure  12)  were  taken  during  the  test  to  measure  the 
time-varying  positions  of  the  tip  vortices  relative  to  the  focal  points. 
Specifically,  they  were  used  to  determine  the  rotor  wake  boundary  in  the 
plane  of  measurements,  and  to  establish  where  the  tip  vortices  crossed 
the  plane  for  specific  blade  azimuth  positions.  The  vortex  positions 
and  their  circulation  sense  for  two  rotor  positions  are  shown  in  Figure 
59*  Coincidentally,  for  the  0.15  advance  ratio  tested,  when  a blade 
was  at  0 = l8o  deg,  the  focal  points  at  0.75R  were  midway  between  the 
blade's  own  tip  vortex  and  the  tip  vortex  at  r/R  = 0.5  from  the  preced- 
ing blade.  When  the  blades  were  at  iff  = 90  and  270  deg,  the  tip  vortex 
from  the  blade  at  90  deg  crossed  the  measurement  plane  at  the  measure- 
ment station,  0.75R.  The  experimental  wake  boundary,  as  defined  by  the 
centers  of  the  tip  vortices  at  their  intersections  with  the  measurement 
plane,  is  shown  in  Figure  58.  Also,  shown  are  the  undistorted  and  analyt- 
ical distorted  wake  boundaries.  Although  slightly  more  contracted,  the 
distorted  wake  boundary,  which  was  calculated  using  the  UTRC  Wake  Geom- 
etry Analysis  described  previously,  agrees  closely  with  the  mean  experi- 
mental wake  boundary  (within  the  accuracy  and  excursions  of  the  test  data), 
Predicted  and  experinental  tip  vortex  distortions  in  the  streamwise 
direction  near  the  measurement  points  were  found  to  be  small  (within 
0.01R)  for  the  selected  test  conditions. 


In  the  selection  of  the  focal  points  for  use  in  the  correlation 
:.tudy,  the  experimental  accuracy  of  the  test  data  was  considered.  Since 
the  LV  test  was  primarily  for  demonstration  purposes,  sophisticated 
instrumentation  to  statistically  process  the  LV  data  was  not  used. 

Taking  advantage  of  the  periodicity  associated  with  the  steady  test  con- 
dition, the  experimental  results  were  obtained  by  fairing  the  data  from 
many  rotor  revolution:;.  In  this  manner,  time  histories  representing  the 
mean  of  the  instantaneous  velocities  were  acquired.  Although  this  pro- 
cedure is  reasonably  accurate  when  the  fluctuations  of  the  data  are 
small,  accuracy  was  lost  at  some  focal  points  due  to  large  fluctuations. 
As  expected,  these  fluctuations  were  generally  the  largest  at  focal 
points  close  to  the  wake  boundary.  The  focal  points  were  thus  selected 
from  those  having  the  least  data  scatter.  It  is  noted,  however,  that 
the  data  fairings  presented  are  believed  to  be  reasonable  but  could 
include  some  inaccuracy,  especially  for  those  points  within  the  wake. 

For  example,  high  harmonic  content  could  possibly  have  been  faired  out. 


6l 


Time  histories  of  the  instantaneous  vertical  velocity  component  at 
three  focal  points,  0.75R  forward  of  the  hub,  are  presented  in  Figure  60. 
The  selected  focal  points  are  0.08R  and  0.02R  above  the  blade  and  0.02R 
below  the  blade,  and  they  represent  points  away  from  the  blade,  near  the 
blade,  and  in  and  out  of  the  wake.  The  velocities  are  presented  for  vary- 
ing time  as  related  to  the  azimuth  position  of  a blade,  and  are  shown  to 
be  periodic  with  blade  passage  interval  (two  cycles  per  rotor  revolution). 
The  instants  of  time  corresponding  to  blade  passages  directly  above  or 
beneath  a focal  point  are  at  & = 0 and  180  deg.  As  described  in  the  ref- 
erence, considering  the  circulation  sense  and  position  of  the  bound  and 
tip  vortex  of  each  blade  relative  to  each  focal  point,  the  flow  direc- 
tions and  locations  of  the  peak  velocities  are  found  to  be  as  expected. 

For  the  points  under  consideration  at  r/R  = 0.75R>  at  a blade  passage 
time,  a selected  focal  point  is  midway  between  the  adjacent  tip  vortices 
as  shown  in  Figure  59*  The  distance  of  the  tip  vortices  and  their  cancel- 
ling effect  allow  the  bound  vorticity  to  predominate.  The  rapid  Lhange  of 
the  velocity  direction  from  upflow  to  downflow,  as  indicated  near  V = 0 
and  180  deg,  substantiates  this.  When  the  blades  are  at  90  and  270  deg, 
a tip  vortex  passes  directly  beneath  the  focal  points.  The  distance  of 
the  blades  and  their  cancelling  bound  vorticity  now  allows  this  tip  vor- 
tex to  predominate,  as  noted  by  the  rapid  change  from  downflow  to  upflow 
before  and  after  the  vortex  passage  (\jr  = 90  and  270  deg  in  the  figures). 
The  combination  of  the  above  two  effects  to  form  a fairly  smooth  time 
history  is  illustrated  in  Figure  6l.  Of  course,  the  relative  degree  of 
influence  of  the  separate  effects  is  dependent  on  the  relative  distances 
from  the  tip  vortex  and  blade  to  the  particular  focal  point.  For  example 
the  more  predominant  influence  of  the  blade  (bound  vortex)  at  the  z'/R 
= 0.08  focal  point  is  indicated  in  Figure  60(a)  by  the  steeper  gradients 
at  .j,  = 0 and  180  deg.  Similarly,  the  influence  of  the  tip  vortex  pas- 
sage is  somewhat  more  predominant  for  the  z'/R  = -0.02  point  between  the 
rotor  and  the  wake  boundary,  as  indicated  in  Figure  60(c)  by  the  steeper 
gradients  at  V = 90  and  270  deg.  As  expected,  the  peak-to-peak  velocity 
amplitudes  increase  with  increasing  proximity  to  the  rotor  and  wake. 

In  Figure  60,  the  time  histories  predicted  using  the  Rotorcraft  Wake 
Analysis  are  generally  in  good  agreement  with  the  time  histories  measured 
with  the  laser  velocimeter.  The  character  and  phasing  of  the  experimental 
time  histories  are  predicted  as  well  as  the  relative  local  gradients  and 
the  growth  in  peak-to-peak  amplitude  with  decreasing  distance  to  the 
rotor  and  wake.  It  is  not  known  whether  the  small  differences  in  peak- 
to-peak  amplitude  and  mean  velocity  are  attributable  to  the  analysis  or 
the  experimental  data,  since  the  differences  are  within  the  accuracy  of 
the  demon strati on  test.  The  increased  peak-to-peak  amplitudes  of  the 
distorted  wake  results  relative  to  those  of  the  undistorted  wake  are  due 
to  the  upward  tip  vortex  distortions  which  place  the  tip  vortex  passage 
closer  to  focal  points  above  it. 
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Time  histories  of  vertical  and  streamwise  instantaneous  velocities 
for  two  focal  points  along  the  vertical  traverse  at  the  r/R  ~ 0.5  location 
are  presented  in  Figures  62  and  63.  The  two  selected  points  consist  of 
a point  above  the  rotor  (at  &z'/r  = 0.04)  and  a point  below  the  rotor  and 
in  the  wake  (at  Az'/r  = -0.04).  These  points  are  shown  in  Figure  58. 

The  character  and  phasing  of  the  vertical  velocities  differ  from  those 
at  the  r/R  = 0.75  location,  mainly  due  to  the  change  in  phasing  of  the 
tip  vortex  passage.  As  shown  in  Figure  59>  the  tip  vortex  passes  these 
points  when  the  blades  are  at  tp  = 0 and  180  deg.  When  the  blades  are  at 
tft  = 90  and  270  deg,  the  blades  and  vortices  are  at  their  furthest  dis- 
tances from  these  focal  points  and  are  in  essentially  self-cancelling 
positions.  Although  the  bound  and  tip  vortex  effects  on  the  vertical 
velocity  are  opposing  before  and  after  their  passage,  they  produce  maximum 
velocity  gradients  near  their  passage  at  ip  = 0 and  180  deg.  The  magnitude 
of  the  streamwise  velocity  agrees  closely  with  the  free-stream  velocity 
(V  = 46  fps)  at  90  and  270  deg,  where,  as  mentioned  above,  the  blades  and 
tip  vortices  have  the  least  influence.  Since  the  bound  vorticity  does  not 
contribute  to  the  streamwise  velocity  component  when  a blade  is  above  or 
below  the  focal  points  under  consideration,  the  increased  streamwise  veloc- 
ities noted  in  Figures  62  and  63  near  i//  = 0 and  180  deg  are  produced  by 
the  tip  vortex  passages.  The  character  and  phasing  of  the  predicted  time 
histories  at  the  0.5R  location  are  generally  similar  to  the  test  data; 
however,  some  differences  do  exist.  In  the  wake  at  Az'/R  = -0.04  (Figure 
3),  the  effect  of  the  tip  vortex  appears  to  predominate  and  produces  down- 
flow prior  to  its  passage  and  upflow  following  its  passage  by  the  points 
at  r/R  = 0.5.  The  closer  passage  of  the  distorted  tip  vortex  produces  a 
higher  predicted  peak  velocity.  The  high  harmonic  content  in  the  predicted 
time  histories  are  believed  to  be  artificially  produced  by  the  finite  fil- 
ament representation  of  the  inboard  wake.  In  Figure  62,  the  velocities 
for  the  point  0.04R  above  the  rotor  are  shown  to  be  reduced"  relative  to 
those  at  the  point  O.OUr  below  (Figure  63),  which  is  consistent  with  the 
increased  distance  from  the  tip  vortices. 


Generally,  it  has  been  shown  that  the  time  histories  of  induced  veloc- 
ities predicted  using  the  Rotorcraft  Wake  Analysis  are  in  favorable  agree- 
ment with  the  limited  amount  of  useful,  time-history  test  data.  Similar 
data  for  further  evaluation  of  the  analysis  at  other  test  conditions  and 
azimuth  positions  are  required. 


Comparison  With  Test  Data  of  Bigger s and  Orloff 


The  use  of  a laser  velocimeter  to  measure  the  flow  velocities  in  the 
wake  of  a two-bladed  model  helicopter  rotor  is  also  reported  in  the 
following  reference  (Reference  26): 


I 
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Diggers,  J.  C. , and  K.  L.  Orloff,  Laser  Velocimeter  Measurements 
of  the  Helicopter  Rotor-Induced  Flow  Field,  Journal  of  the  Ameri- 
can Helicopter  Society,  Vol.  20,  No.  1,  January  1975* 


Laser  velocimetry  techniques  were  used  to  measure  two  components  of  induced 
velocity  in  the  outboard  vicinity  of  blades  at  two  azimuth  positions. 
Time-averaged  and  instantaneous  velocity  components  in  the  vertical  and 
streamwise  directions  (vz<  and  vx«)  were  measured  in  the  NASA-Ames  7- 
by  10-foot  wind  tunnel  at  a rotor  advance  ratio  (v/f)R)  of  0.18.  Measure- 
ments were  made  at  a series  of  stationary  poixits  in  a vertical  plane 
through  the  center  of  the  hub  (the  x-y  plane)  at  the  advancing  and 
retreating  blades  0 = 90  and  270  deg)u  Results  were  obtained  primarily 
for  the  advancing  blade.  Time- averaged  velocities  were  measured  by  con- 
tinuously averaging  the  velocities  of  seeded  particles  passing  through 
the  fixed  focal  points.  To  measure  instantaneous  velocities,  a strobe 
technique  was  used  which  relies  upon  the  periodicity  of  the  rotor  flow 
for  steady-state  conditions.  Data  were  accepted  only  once  per  revolution 
during  a short  "data  window",  which  was  open  only  when  the  blade  was  at  a 
prescribed  azimuth.  Information  on  the  rotor  design  and  single  test  condi- 
tion is  presented  in  Table  11.  Thrust  was  not  measured. 


TABLE  11.  ROTOR  AND  TEST  INFORMATION  --  BIGGERS  & ORLOFF 


Single,  model  , hingeless  rotor 


Rotor  Desit 


Number  of  blades 


Blade  taper  ratio 


Rotor  radius,  ft 
Blade  chord,  ft 
Rotor  solidity 


0.354 


Blade  aspect  ratio 
Blade  twist,  deg 


0.0644  Airfoil  section 


NACA  0012 


Test  Condition 

/i.  = 0.18,  OR  = 110  fps,  9 ^ = 8.5°,  ct  = -6.6°,  vq  = 1.83  fps 


A^  = 03,  B^  = (?,  not  measured  (calculated  = 0.0048) 
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In  order  to  show  the  sensitivity  of  the  predicted  velocities  to  tip 
vortex  distortions,  the  Rotorcraft  Wake  Analysis  was  run  with  both 
undistorted  and  distorted  wake  representations.  The  test  control  angles 
and  rotor  tilt  angles  were  used  in  the  analysis.  Top,  side,  and  rear 
views  of  the  analytical  undistorted  and  distorted  tip  vortex  geometries 
relative  to  the  tip-path  coordinate  axes  are  presented  in  Figures  64  and 
65,  respectively.  The  undistorted  wake  skew  angle,  X,  is  78.2  deg.  The 
tip  vortices  from  both  blades  are  shown  for  the  rotor  position  at  which 
instantaneous  velocities  were  measured  at  the  advancing  and  retreating 
blades.  By  comparing  the  tip  vortex  positions  between  these  two  figures, 
it  is  evident  that  the  tip  vortex  distortions  predicted  by  the  UTRC  Wake 
Geometry  Analysis  (Deck  F506)  are  relatively  small  near  the  blades  for 
this  rotor  position  and  test  condition.  Greater  tip  vortex  distortions 
occur  downstream,  but  they  are  of  minor  consequence  to  the  induced  veloc- 
ities near  the  advancing  and  retreating  blades. 

The  location  of  the  tip  vortex  from  the  preceding  blade  was  deter- 
mined during  the  test  when  it  crossed  the  advancing  blade  at  \fj  = 90  deg. 

A comparison  of  the  test  tip  vortex  location  with  the  analytical  undis- 
torted  and  distorted  tip  vortex  locations  is  presented  below; 


Measured  tip  vortex  location 
Undistorted  tip  vortex  location 
Distorted  tip  vortex  location 


x'/R 

y’/R 

z’/R 

Az  '/R 

0.0 

0.85 

-0.05 

-0.048 

0.0 

0.87 

-0.104 

-0.102 

0.0 

0.86 

-0.045 

-0.043 

The  z'/R  values  are  relative  to  the  rotor  hub  center  in  the  wind  tunnel 
coordinate  system.  Considering  the  preconing  and  lateral  tilt  of  the 
rotor,  the  vertical  distances  between  the  tip  vortex  and  the  blade  are 
listed  as  Az'/R.  As  indicated  above,  the  predicted  and  measured  tip 
vortex  positions  are  in  good  agreement. 

The  velocities  calculated  using  the  Rotorcraft  Wake  Analysis  with 
undistorted  and  distorted  wake  representations  are  compared  with  the 
laser  velocimeter  test  data  in  Figures  66  through  71.  In  most  cases, 
the  plot  format  of  the  reference  has  been  retained  to  maintain  a direct 
relationship  with  the  test  data  as  presented  therein.  Velocities  are 
nondimen sionali zed  by  tip  speed,  and  the  wind  tunnel  axis  system  is  used. 
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Time-Averaged  Velocities 

In  Figure  66,  the  predicted  and  experimental  radial  distributions  of 
the  time-averaged  vertical  velocities  are  presented  for  several  vertical 
positions  near  the  rotor.  The  vz»  axis  for  each  plot  is  placed  at  the 
z'/R  location  for  that  set  of  data  (the  z'/R  scale  is  expanded  to  ten 
times  the  y'/R  scale  for  clarity  of  presentation).  The  radial  distribu- 
tions are  characterized  by  downflow  inside  the  wake,  changing  to  upflow 
outside  the  wake  due  to  the  passage  of  the  tip  vortices.  The  gradients 
diminish  with  increasing  distance  from  the  rotor.  The  analytical  velocity 
distributions  are  in  agreement  with  these  features.  The  variations  of 
time-averaged  vertical  velocities  (predominantly  downwash)  with  vertical 
position  are  plotted  in  Figure  6?  for  several  radial  locations.  The  vz> 
axis  for  each  traverse  is  placed  at  the  y'/R  location  for  that  data  set. 
Again,  the  predicted  velocity  distributions  exhibit  the  same  character  as 
the  measured  values,  although  absolute  differences  exist.  Moving  down 
each  vertical  traverse,  the  predicted  and  measured  downwash  increases  to 
a peak,  decreases  rapidly  to  the  vicinity  of  the  time-averaged  wake 
boundary,  and  then  gradually  diminMies  to  zero.  Differences  in  location 
of  the  primary  velocity  gradients  in  Figures  66  and  67  indicate  the 
extent  of  the  difference  between  the  predicted  and  measured  time-averaged 
wake  boundaries.  In  Figure  68,  the  predicted  and  measured  radial  distri- 
butions of  the  time-averaged  streaniwise  velocity  are  compared  at  the 
z'/R  = -0.05  location  (the  free-stream  V/OR  = 0.18  is  included). 

Instantaneous  Velocities 


Radial  distributions  of  instantaneous  vertical  velocities  are  pre- 
sented in  Figure  69  for  the  same  vertical  positions  as  Figure  66.  As 
mentioned  above,  the  measured  tip  vortex  location  was  at  x'/r  = O.85, 
z'/R  = -0.05.  This  tip  vortex  location  is  evident  in  the  measured  instan- 
taneous velocity  profile,  which  shows  a sharp  downward  peak  at 
y'/R  = 0.845  and  a less  pronounced  upward  peak  at  y'/R  = 0.860.  The  pre- 
dicted instantaneous  velocities  based  on  the  di storted  wake  representation 
are  much  improved  over  the  undistorted  wake  results  due  to  the  more  accu- 
rate positioning  of  the  tip  vortex.  At  the  z'/R  = -0.05  elevation,  the 
presence  of  the  distorted  tip  vortex  is  also  evident  (exact  location  is 
z'/R  = -0.045),  whereas  the  undistorted  tip  vortex  is  most  evident  at 
z'/R  = -0.09  (exact  location  is  z'/R  = -0.104).  Consistent  with  the 
experimental  results,  the  influence  of  the  predicted  distorted  tip  vortex 
from  the  preceding  blade  is  less  pronounced  as  the  distance  away  from 
the  core  increases.  However,  its  influence  is  still  significant  as  far 
as  z'/R  = -0 .09. 
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A comparison  of  the  radial  distributions  of  instantaneous  stream- 
wise  velocity  for  the  z'/R  = -0.05  elevation  is  included  in  Figure  68. 

The  limitations  of  the  Rotorcraft  Wake  Analysis  for  predicting  the  local 
instantaneous  velocities  in  the  immediate  vicinity  of  blades  with  a high 
degree  of  accuracy  are  illustrated  in  Figures  70  and  71.*  It  will  be 
shown  that  the  limitations  are  primarily  due  to  the  lifting  line  repre- 
sentation of  the  blades  in  the  analysis.  However,  before  describing  the 
results,  it  is  noted  that  when  the  analysis  is  uned  to  compute  induced 
velocities  at  the  blades  for  blade  airloads  and  performance  purposes, 
the  influence  of  the  blade  bound  vortex  (lifting  line)  at  blade  segment 
points  along  itself  is  eliminated.  Although  it  is  recognized  that  finite 
chord  and  thickness  effects  on  the  local  flow  at  the  blade  are  neglected, 
the  oversensitivity  of  the  local  flow  predictions  to  small  distances  from 
a concentrated  oound  vortex,  such  as  will  be  shown,  does  not  exist  in  such 
applications. 

In  Figure  70,  the  predicted  and  measured  instantaneous  vertical, 
velocities  at  r/R  = 0.95  are  presented  for  the  advancing  and  retreating 
blades,  respectively.  In  the  reference,  the  velocity  profile  of  the 
experimental  data  in  Figure  70  (a),  showing  upward  flow  above  the  blade  and 
downward  flow  below  it,  is  attributed  primarily  to  the  influence  of  the 
airfoil  and  its  attitude;  that  is,  (l)  streamline  divergence  effects  due 
to  airfoil  profile  (maximum  thickness  is  at  the  30-percent  chord,  and  the 
measurements  were  made  at  the  25-percent  chord),  and  (2)  the  upward  com- 
ponent of  blade  circulation  at  the  upper  surface,  and  the  downward 
component  at  the  lower  surface  associated  with  the  1.5-deg  nosedown 
attitudes  of  the  airfoil  section  relative  to  the  measurement  axes. 

These  effects  are  neglected  in  the  analysis.  In  fact,  the  bound  vortex 
in  the  analysis  does  not  affect  the  vertical  velocities  at  all  for  these 
results  because  the  vertical  axes  in  Figure  70  pass  directly  through  the 
quarter-chord  where  the  bound  vortex  is  located.  The  analytical  velocity 
profile  in  Figure  70(a)  is  attributed  to  the  influence  of  the  tip  vortices. 
That  is,  as  the  traverse  is  made  starting  above  the  blade,  increasing  down- 
flaw  is  due  to  the  advancing  blade's  own  tip  vortex  and  then  the  upflow 
is  due  to  the  tip  vortex  of  the  preceding  blade  which  is  at  a lower  posi- 
tion. For  a similar  traverse  for  the  retreating  blade  shown  in  Figure 
70(b),  the  velocity  profile  is  different.  The  peak  downflow  at  the  blade 
may  be  attributed  to  the  retreating  blade's  own  tip  vortex,  and  the  lack 
of  a similar  upflow  peak  may  be  attributed  to  the  position  of  the  previous 
blade's  tip  vortex  at  a much  larger  distance  from  0.95R  on  the  retreating 
side  (see  Figure  65).  However,  the  absence  of  experimental  upflow  below 


*Vertical  distances  in  Figures  70  and  71  are  relative  to  the  blade 
(Az’/lO  rather  than  the  hub  (z/r)  as  ir.  previous  figures. 
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the  blade  due  to  the  above  two  real  airfoil  influences,  cited  in  the 
reference,  is  somewhat  contradictory,  particularly  since  the  airfoil  is 
at  a higher  noseup  attitude  (l8  deg)  relative  to  the  x'  axis. 


Vertical  traverses  of  the  instantaneous  streamwise  velocities  near 
the  advancing  and  retreating  blades  at  r/R  = 0.95  are  presented  in 
Figure  71.  . Since  the  blades'  own  tip  vortices  are  essentially  aligned 
with  the  x'  axis  near  the  measurement  plane  when  the  blades  are  at  the 
azimuth  positions  investigated,  their  contribution  to  the  streamwise 
velocities  is  small.  Thus,  the  bound  vorticity  is  the  major  contribution 
to  the  streamwise  velocities  in  these  figures.  In  both  figures  the  pre- 
dicted velocities  result  mainly  from  the  bound  vortex  induced  flow  toward 
the  trailing  edge  above  the  blade,  and  toward  the  leading  edge  below  the 
blade.  This  induced  flow  is  superposed  on  the  free- stream  flow.  For  the 
retreating  blade  (Figure  71(b)),  the  experimental  flow  reversal  is  evident. 
However,  for  the  advancing  blade  (Figure  71(a)),  it  is  not.  In  the 
latter  figure,  the  analytical  induced  velocities  are  accurate  above  the 
blade.  Just  below  the  blade,  the  reason  for  the  absence  of  upstream 
induced  velocities  (relative  to  the  free-stream  velocity)  in  the  test 
data,  expected  from  bound  vorticity,  is  not  known.  The  magnitudes  of  the 
predicted  velocity  peaks  in  Figure  71(b)  are  dependent  on  the  lifting 
line  representation,  and  they  were  found  to  be  very  sensitive  to  small 
changes  in  chordwise  variations  of  the  traverse  line.  The  peak  values 
and  their  location  are  also  sensitive  to  the  artificial  vortex  core  size 
(diameter  = 0.15  chord)  assigned  to  the  bound  vortex  to  preclude  veloc- 
ities approaching  infinite  magnitude  very  close  to  the  lifting  line.  In 
addition,  the  accuracy  of  the  velocity  values  at  a point  near  the  blade 
is  dependent  on  the  accuracy  of  the  calculated  circulation  strengths  of 
the  bound  vortex  segments  near  the  point. 

It  is  noted  in  Figures  70  and  71  that  the  distorted  wake  results  are 
similar  to  the  un distorted  wake  results.  This  is  explained  by  the  rela- 
tively small  distortions  and  the  small  effect  of  the  tip  vortices  of  the 
preceding  blades  on  the  vertical  and  streamwise  instantaneous  velocities 
at  the  0.95R  locations  for  these  figures.  Greater  wake  distortions  and 
associated  induced  effects  occur  at  other  blade  azimuths  than  those  tested. 

Generally,  the  results  presented  above  indicate  that  the  Rotorcraft 
Wake  Analysis  predicts  the  features  of  the  velocity  distributions  in  the 
wake,  but  the  degree  of  accuracy  varies  depending  on  the  accuracy  of  the 
wake  representation.  Also,  as  expected,  for  points  removed  from,  but  in 
the  immediate  vicinity  of,  a passing  blade,  a lifting  line  representation 
of  the  blade  is  insufficient.  Laser  velocimetry  is  shown  to  be  an  excel- 
lent means  of  providing  data  for  correlation  purposes.  LV  data  for  other 
test  conditions  and  at  other  blade  azimuths  where  wake  dirtortions  are  more 
significant  are  required. 
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Comparison  With  Test  Data  of  Heyson,  KACA  TN  3242 


Following  the  incorporation  in  the  Rotorcraft  Wake  Analysis  of  a 
provision  for  the  computation  of  off-rotor  velocities  for  dual  rotor.  , 
the  computer  analysis  was  applied  to  calculate  the  induced  velocity 
field  of  the  tandem-  and  single-rotor  configurations  described  in  the 
following  reference  (Reference  27): 

e/ron.  !:.  i'rel  imi  :.ary  ! ee  ;1  ' s from  Low-Field  "easurements 

Around  .'i:.sLe  a:.d  i'andem  i-.o*  ors  is  'he  Lan.-le.y  Full  -:.‘cale 
el  . FA  "A  7747.  er  Is  *54. 

The  rotor  and  test  information  from  this  reference  are  presented  in  Table 
12.  The  15-foot  diameter  rotor,  were  tested  in  the  NAOA-Langley  Full - 
.'cale  Wind  Tunnel  at  an  advance  ratio  of  0.15.  One  of  these  rotors  was 
used  in  the  previously  described  experimental  program  of  Reference  24. 

A five-tube  pressure  survey  rake  was  used  to  measure  local  downwash 
angles  in  the  :urvey  planes  indicated  in  Figure  77.  The  local  downwash 
angles  were  averaged  spanwi.e  fi.e.,  laterally)  at  each  elevation  and 
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nondimensionalized  by  ihe  downwash  angle  at  the  center  of  a single-rotor 
based  on  momentum  theory.  Single-rotor  data  were  acquired  by  simply 
removing  the  rear  rotor  from  the  tandem  configuration.  The  tandem  rotor 
was  set  to  operate  at  the  same  lift  coefficient  as  the  single  rotor  by 
maintaining  the  single-rotor  controls  for  the  forward  rotor  and  adjusting 
the  rear  rotor  controls.  These  controls  were  not  provided  in  the  refer- 
ence, but  were  estimated  to  correspond  to  the  specified  lift  coefficient  I 

c 

and  wake  momentum  skew  angle  (82.7°).  ^ 

i 

Comparisons  of  the  predicted  results  from  the  Rotorcraft  Wake  Anal-  3 

ysis  with  the  experimental  data  are  presented  in  Figure  72.  As  noted,  | 

the  trends  with  vertical  and  longitudinal  position  and  the  trends  from 
single  to  tandem  rotors  are  accurately  predicted.  Good  agreement  with 


test  data  is  shown  for  both  the  single  and  tandem  rotors  in  the  forward 
survey  plane.  However,  the  accuracy  deteriorates  in  the  rearward  survey 
planes.  This  inaccuracy  is  believed  to  be  inherent  for  the  undistorted 
wake  model  used  in  the  analysis.  The  wake  downstream  from  a rotor  exhib- 
its significant  distortions  which,  as  shown  in  the  reference,  are  charac- 
terized by  the  downstream  transformation  of  the  wake  into  a single  pair 
of  vortices  for  both  single  and  tandem  rotors.  Thus,  significant  dis- 
tortions of  the  wake  from  the  forward  rotor  mixed  with  the  wake  from  the 
rear  rotor  occur  near  the  rear  rotor,  "’hese  are  not  represented  in  the 
wake  model.  Also,  some  inaccuracy  may  be  associated  with  the  uncertain- 
ties in  the  controls  and  dependency  of  the  spanwise  average  on  the 
selection  of  points,  particularly  near  the  wake  boundaries.  Local 
velocity  data  would  have  been  more  desirable  for  comparison  purposes. 
Nonetheless,  these  initial  results  indicate  that  improvements  in  the 
tandem-rotor  downwash  predictions  may  be  anticipated  through  the  use 
of  the  dual  rotor  provision  in  the  Rotorcraft  Wake  Analysis,  which 
accounts  for  the  influence  of  the  wakes  of  both  rotors. 
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CONCLUDING  REMARKS 

The  development  of  the  UTRC  Rotorcraft  Prescribed  Wake  Induced 
Velocity  Analysis  ( "Rotorcppft  Wake  Analysis")  provides,  in  a single  com- 
puter analysis,  the  capability  for  calculating  flow  velocities  induced 
by  the  rotor (r, ) for  the  following  field  points,  configurations,  condi- 
tions, and  wake  representations: 

a.  Time -aver aged  and  instantaneous  induced  velocities  (three 
components)  on  and  off  the  rotor. 

b.  Single  and  dual  rotor  configurations  (coaxial,  tandem,  side- 
by-side,  etc.). 

c.  Steady  hover  and  forward  flight  conditions. 

d.  Undistortcd  and  distorted,  analytical  and  experimental,  pre- 
scribed wake  geometries. 


For  rotor  configurations  and  conditions  for  which  test  data  were 
available,  the  induced  velocities  calculated  in  the  rotorcraft  wake 
analysis  are  generally  in  good  agreement  with  measured  velocities.  The 
successful  prediction  of  primary  characteristic  features  of  both  time- 
averaged  and  instantaneous  induced  velocity  distributions  in  and  out  of 
the  rotor  wake,  and  local  velocity  variations  with  several  rotor  design 
and  test  condition  parameters,  has  been  demonstrated  within  the  scope 
of  the  available  test  data.  However,  the  degree  of  accuracy  at  a speci- 
fic field  point  is  sensitive  to  the  accuracy  of  the  prescribed  wake  rep- 
resentation, proximity  of  the  field  point  to  a wake  boundary  or  blade, 
and  the  input  flight  condition,  rotor  control,  and  blade  motion  informa- 
tion. 
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RECOMMENDATIONS 


The  Rotorcraft  Wake  Analysis  represents  the  combination  and  exten- 
sion of  component  computer  programs  which  have  been  developed  and  in  use 
at  the  United  Technologies  Research  Center  for  many  years.  During  this 
time,  the  component  analyses  have  been  applied  to  many  and  varied  prob- 
lems, and  the  results  have  generally  been  favorable.  However,  like  any 
analysis,  the  Rotorcraft  Wake  Analysis  should  continue  to  be  refined 
and  checked  to  remove  limitations  and  to  incorporate  advances  in  appli- 
cable technology  as  they  become  available.  The  following  recommenda- 
tions for  acquiring  additional  experimental  data,  performing  a sensitiv- 
ity study,  continuing  application  and  evaluation,  and  providing  refine- 
ments to  the  analysis  are  thus  offered: 

1.  Additional  Experimental  Data  - Additional  experimental 
data  are  required  to  further  evaluate  the  Rotorcraft  Wake  Analysis, 
and  to  provide  wake  geometry  data  for  improving  the  analytical  rep- 
resentation. Systematic  data  should  be  acquired  over  an  increased 
range  of  rotor  designs,  flight  conditions,  and  field  points  for 
single-rotor,  dual-rotor,  and  winged  configurations  of  interest.  In 
addition  to  flow  velocity  data,  wake  geometry  and  airloads  data 
should  be  acquired  concurrently  to  provide  a consistent  set  of 
correlative  information.  Laser  velocimeter  techniques  should  con- 
tinue to  be  used  to  provide  velocity  data  at  the  rotor  blades  as 
well  as  in  the  wake,  and  to  provide  blade  airload  distributions 
(bound  vorticity) . The  analysis  should  be  used  for  guidance  in  a 
systematic  selection  of  test  conditions  and  field  points  to  acquire 
the  most  useful  data  cost-effectively. 
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2.  Sensitivity  Study  - The  sensitivity  of  the  induced  veloc- 
ities predicted  by  the  Rotorcraft  Wake  Analysis  to  small  variations 
in  input  flight  condition,  rotor  controls,  blade  motions,  and  wake 
representation  should  be  investigated  for  selected  sources  of  cor- 
relative data  to  determine  the  degree  of  accuracy  required  of  the 
input  parameters  and  test  information  used  to  establish  the  input. 
The  sensitivity  of  the  correlation  results  to  the  assumed  wake  rep- 
resentation should  be  investigated  by  systematically  varying  the 
wake  parameters  (wake  azimuth  increment,  number  of  wake  revolutions, 
vortex  core  size,  vortex  filament  distribution,  etc.)  in  the  analy- 
sis. This  has  been  accomplished  to  seme  extent  in  the  past  with 
the  component  programs  of  the  Rotorcraft  Wake  Analysis  by  comparing 
rotor  performance  and  blade  bending  moments.  With  the  availability 
of  additional  test  data  at  and  away  from  the  rotor  and  the  provision 
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of  the  off-rotor  velocity  capability  in  the  analysis,  a more 
accurate  sensitivity  study  should  be  performed  to  assist  with 
future  applications  and  evaluations  of  the  analysis. 

3.  Continued  Application  and  Evaluation  - Application  of  the 
Kotorcraft  Wake  Analysis  in  conjunction  with  rotor  performance  and 
blade  response  analyses,  and  comparison  of  the  results  with  test 
data,  should  be  continued  to  further  evaluate  the  induced  veloci- 
ties, airloads,  rotor  performance,  blade  response,  and  bending 
moments  predicted  with  the  combined  analysis.  As  correlative  test 
data  become  available,  the  range  of  configurations  and  test  condi- 
tions for  which  induced  velocities  are  compared  should  be  extended. 
In  particular,  further  comparison  is  required  for  low-speed  test 
conditions  and  dual-rotor  and  winged  helicopter  configurations. 

4.  Refinements  ~ The  Rotorcraft  Wake  Analysis  has  been 
organized  to  allow  components  applicable  to  specific  technical  areas, 
configurations,  or  flight  mode  to  be  updated  or  replaced  as  the 
state  of  the  art  advances  or  tie  need  arises.  Also,  the  use  of  a 
prescribed  wake  in  the  analysis  retains  flexibility  in  the  type 

and  source  of  wake  geometry  representation.  As  additional  accurate 
and  generalized  wake  geometry  information  becomes  available  from 
experimental  programs  and/or  wake  geometry  analyses,  it  can  be  in- 
corporated in  the  analysis.  Wake  modeling  refinements  to  account 
for  wake  mixing,  vortex  core  dissipation,  and  inboard  wake  distor- 
tions should  eventually  be  incorporated.  The  inclusion  of  an 
optional  replacement  of  the  lifting  line  representation  of  the  blade 
by  a lifting  surface  representation,  at  least  in  local  areas  of 
significance,  should  be  considered.  Also,  the  analysis  should  be 
extended  to  include  a fuselage  and  tail  rotor,  operation  near  the 
ground  (in  ground  effect),  and  operation  in  transient  (maneuvers 
and  unsteady  wind)  condition.*- . 
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Figure  15.  Flow  Chart  of  the  UTRC  Prescribed  Wake  Hover  Perfor- 
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Figure  20.  Parameters  for  Wing  Induced  Aerodynamic  Inter- 
ference at  a Rotor  . 
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Figure  22.  Analytical  Normal  Velocity  Component  Induced  by  Wing 
at  Rotor  . 
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Figure  24.  Up  Vortex  Positions  for  the  Three  Test  Conditions 
of  Boatwright,  USAAMRDL  TR  72-33,  n - 0. 
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Figure  26  Comparison  of  Instantaneous  Normal  Component  of 

Induced  Velocity  at  Several  Distances  from  Rotor — 
Boatwright,  USAAMRDL  72-33;  u =»  0,  Test  Condition 
1,  ■ 0 deg. 
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Figure  28.  Comparison  of  Radial  Distributions  of  Instantaneous 

Radial  Velocity  Component  for  Three  Rotor  Test  Condi- 
tions and  Several  Distances  Under  Rotor — Boatwright, 
USAAMRDL  TR  72-33,  n*  0,  >J>  » 0 deg. 
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Figure  32.  Comparison  of  Time-Averaged  Nonnal  Component  of 
induced  Velocity  for  Two  Tip  Speeds  and  Classical 
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z/R  a -0.145  and  -0.45.  ^ ’ 


ROTOR  DISK, 


RADIAL  COORDINATE,  r/R 


Figure  37.  Experimental  Hme-Averaged  Normal  Induced  Velocity  Con- 
tours --  Fradenburgh  Model  Rotor,  fi  = 0 . 
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Figure  38.  Analytical  Time-Averaged  Normal  Induced  Velocity  Con 
tours  Based  on  Classical  Wake  Representation — Fraden 
burgh  Model  Rotor,  0 = 0 . 
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Figure  49.  Comparison  of  Time-Averaged  Normal  Component  of  Induced 
Velocity  Along  the  Centerlines  of  the  Rotor  Disk  — 
Heyson  and  Katzcff,  NACA  TR  1319)  Condition  2, 
g a 0.l40,  z/R  a 0 . 
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Figure  50.  Comparison  of  Time-Averaged  Normal  Component  of 
Induced  Velocity  for  x/R  - 0 and  0. 5 — Heyson 
and  Katzoff,  NACA  TR  1319,  Condition  1,  „ . 0.095. 


VELOCITY  COMPONENT  NORMAL  TO  ROTOR 


VELOCITY  COMPONENT  NORMAL  TO  ROTOR  DISK.  -v-/v 


(WITH  SMALL  SYMBOLS!  MEASURED  INACA  TR1319I 

ROTORCRAFT  WAKE  ANALYSIS 
(UNDISTORTED  WAKE) 


PRESSURE 
'/R  TESJ  DATA 

0.22  O 

-0.26  □ 

0.49  A 


HOT  FILM 
TEST  DATA 

d 

d 

a 


ROTORCRAFT  WAKE  ANALYSIS 
_ (CL  ASSICAL  WAKE) 


RADIAL  COORDINATE.  r/R 


Figure  54.  Comparison  of  Radial  Distributions  of  lime-Averagec 
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Figure  60.  Comparison  of  Time  Histories  of  the  Vertical  Com- 
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Figure  62.  Comparison  of  the  lime  Histories  of  Vertical  and 

Streamwise  Velocity  for  a Focal  Baint  at  r/p,  ® 0.50 
dz'/R  » 0.04  — Landgrebe  and  Johnson,  n - 0.15 . 
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Figure  67.  Comparisons  of  Vertical  Distributions  of  Hme- 
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APPENDIX  A 

BIOT- SAVANT  LAW  AS  APPLIED  TO  A VORTEX  ELEMENT 


The  velocities  induced  at  an  arbitrary  point  P by  a straight  vortex 
element  of  strength  r and  bounded  by  end  points  A and  B can  be  computed 
using  the  classical  Biot-Savart  Law.  The  equations  presented  below  are 
in  Cartesian  form.  These  equations  apply  only  in  the  irrotational  flow 
region  outside  of  the  vortex  core. 


= (^r)s  [(yp’yA)(Zp-zB)  - (zP-2A)(yP-yB) 

| (A-l) 

(47r)^  (^Zp  ZA^XP~Xb)_  ^ XP  ~ XA^ZP  ~ZB) 

1 (A-2) 

: |(xp  'xA)(yP  -ya)-  (yP-yA)  (xp"xb)] 

(A-3) 

- 1 (AP  + BP)/  (AP)  (BP) 

R (AP)(BP)  + I + J + K 

(A-U) 

I - (Xp-XA)  (Xp  -Xg) 

(A-5) 

*■'. = (yp*yA)(yP  -yB) 

(A-6) 

* " (Zp  _Za)  (ZP  ~zB) 

(A-7) 

yj p “ xAf+(yp  -yA)  + (zp  ~zA)Z 

(A-8) 

BP  = ~JWP  -xB)z+(yp  -yB)*+(zP  -zB)2 

(A-9) 
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APPENDIX  B 


PRIMARY  EQUATIONS  FOR  THE  UTRC 
PRESCRIBED  WAKE  hOTOR  INFLOW  ANALYSIS* 

In  order  to  relate  to  the  program  input-output  symbols  and  to  des- 
cribe some  of  the  symbols  generated  within  the  program,  program  symbol 
notation  is  generally  used  in  the  following  equations. 

The  blade  circulation  at  any  blade  segment  I,  for  any  azimuthal 
position  J,  is,  from  the  Kutta-Joukowski  relation, 


where 


CIRC(I,J)  = C(I)  • U(I,J) 


CL(l,J)/2 


CIRC (I, J) 
C(I) 
U(I,J) 
CL(I,J) 


= local  blade  element 
= local  blade  element 
= local  blade  element 
- local  blade  element 


circulation,  fu  /sec 
chord,  ft 

resultant  velocity,  fps 
lift  coefficient 


(B-l) 


The  local  blade  element  resultant  velocity  U(I,J)  can  be  expressed  as 


U(I,J)  = [UT(I,J) +UTl(l,j)r  + [Wl(l,j)  + W(I,J)]C 


(B-2) 


where  UT(I,J)  = OMGR  * ( RS ( I ) + MU*SIN(PSl)) 


(B-3) 


and 


OMGR  = 
MU  - 
RS(I)  = 

UTl(l, j)  = 
W(I,J)  = 


Rotor  tip  speed,  fps 
Rotor  advance  ratio 

Rotor  blade  segment  radial  position  ( nondimens ional- 
ized  by  R) 

Induced  tangential  velocity,  fps 
Induced  inflow  velocity,  fps 


and  with  small  angle  assumptions, 


* The  fundamental  background  for  the  analytical  formulation  used  in 
the  computer  program  was  obtained  from  the  technical  approach  described 
by  Piziali  and  DuWaldt  in  Reference  l4.  Although  a new  computer  program 
has  been  developed,  and  many  modifications  and  refinements  have  been  in- 
corporated, practicularly  in  regard  to  wake  representation,  many  of  the 
fundamental  assumptions  are  retained  frcm  that  reference. 


i 


\ 


155 


Wl(l,j)  = Noninduced  inflow  velocity  (flapping,  pitching, 
aeroelastics , etc,),  fps. 

For  helicopter  rotors,  UTl(l,J)  may  be  neglected.  Also  since  W(l,J)  and 
Wl(l,j)  are  small  in  the  important  blade  tip  region  relative  to  the  ro- 
tational velocity,  it  is  assumed  that 

U(l , J)  = UT(I,J)  (B-U) 

so  that  a linear  closed  form  solution  can  be  obtained.  Also  in  order  to 
linearize  the  solution,  the  lift  coefficient  is  expressed  as 

CL(I,J)  = AT( I , J ) * ALPHAE(I,J)  (B-5) 

where  AT(l,J)  = local  blade  element  lift  curve  slope 

ALPHAE(l,J)  = local  blade  element  angle  of  attack 

In  stall,  the  maxi  mum  lift  coefficient  is  limited  in  the  program  to 
CU4AX(l,J)  which  is  determined  from  input  values  of  angle -of- at  tack 
values  (ALMAX(l,J)  beyond  which  stall  is  assumed.  The  angle  of  attack 
is,  for  small  angle  assumptions  and  consistent  with  the  above  assumptions, 

ALPHAE ( I , J ) = [W1(I,J)  + W(I,J)]  /UT(I,J)  (B-S) 

Substituting  Equations  (B-U)  through  ( B— 6)  into  Equation  (B-l),  the  cir- 
culation equation  becomes, 

CIRC(l, J)  = C(I)  • AT(I,J)  • rwi(l,j)  + W(l,j) ] /2  (B-7) 

Since  there  are  MSIZE  unknowns , 

where  MSIZE  - (ITOT  + 1)  * JTOT  and  ITOT  and  JTOT  are  the  total- 

number  of  blade  segments  and  azimuth  positions, 
there  are  MSIZE  equations  of  the  form  of  Equation  (B-7).  This  set 
of  MSIZE  equations  can  be  written  in  matrix  notation  as: 

CIRC(l,  J)  = rC(l)  • AT(  I,  J ) • (W1(I,-J)  + W(I,J)}]  (B-8) 

The  axial  induced  velocity  in  Equations  (B-7)  and  (B-8)  is  ex- 
pressed, using  the  Bi:t-Savart  law,  as  a function  of  the  trailing  wake 
filaments  and  circulation  strengths.  The  axial  velocity  induced  by  a 
given  wake  segment  (subscript  M)  of  filament  K at  a given  blade  segment 
(l,j)  is  expressed  in  terms  of  the  Biot-Savart  law  as 
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W(K,M)  = CIRCT(K,M)/4ttR  * GC(K,M) 


(B-9) 


CIRCT(K,M)  = 

gc(k,m)  = 


circulation  strength  of  the  Mth  segment  of 
the  Kth  filament,  ft2/sec 

geometric  influence  coefficient  from  the  Biot- 
Savart  law 


The  geometric  coefficient,  GC,  is  a function  of  the  prescribed  wake 
coordinates  X(K,M),  Y(K,M),  Z(K,M).  The  lengthy  expressions  for  GC  may 
be  obtained  from  the  program  listing.  Since  the  circulation  distribu- 
tion is  assumed  to  be  periodic  in  time  for  the  forward  flight  condition, 
each  wake  filament  has  a periodic  circulation  strength  distribution  along 
its  length.  The  total  axial  induced  velocity  due  to  contributions  of 
all  trailing  circulations  for  all  of  the  filaments  at  the  blade  point 
(l,J)  is  the  sum  of  Equation  (B-9)  evaluated  for  all  segments  in  the 
wake.  The  resultant  expression  frcrn  Equation  (B-9)  becanes 


KLIM  MLIM 

W(I,J)  = - E E CIRCT(K,M) 


GC(K,M)  ( B-10 ) 


where  KLIM  is  the  total  number  of  trailing  filaments,  and  MLIM  is  the 
total  number  of  segment  end  points  for  each  filament. 


However  since  the  circulations  are  periodic  in  M this  relationship  can 
be  expressed  by  combining  elements  with  like  circulations, 


W(!,J)  = n- 


KLIM  JTOT 
E E 
K MM 


circt(k,mm)  . gctrow(k,mm) 


(B-ll) 


where  MLIM 

GCTROW(K,MM)  = E GC(K,M) 

(M=MM) 

MM  combines  all  elements  of  a vortex  filament  shed  fr  a blade 
at  the  same  blade  azimuth  position,  and  thus  combines  elements  with  like 
circulations.  The  total  induced  velocity  due  to  all  bound  circulations 
of  the  rotor  can  be  expressed  as, 


WB(I,J)  = 


IT'w  1’  JioT 
E E 
I MM 


;(I,MM)  • GCB(I,MM)  (B-12) 


where  GCB(l,MM)  = geometric  coefficient  due  to  bound  circulation 

(using  Biot-Savart  law) 

WB(l,J)  = induced  axial  velocity' due  to  bound  circulation 
IB  = number  of  blades 
ITOT  = number  of  bound  segments  per  blade 

In  order  to  express  the  total  axial  induced  velocity  due  to  all  of  the 
trailing  circulations  in  terms  of  the  bound  circulation,  the  circulation 
strength  of  each  trailing  filament  segment  may  be  expressed  as  the  dif- 
ference between  the  bound  circulation  of  the  two  blade  segments  adja- 
cent to  its  point  of  origin  at  the  blade  at  the  time  of  generation. 

Thus , 


circt(k,mm)  = CIRC(K,MM)  - CIRC(K-1,MM)  (B-13 ) 

Substituting  Equation  (b-12)  and  Equation  (B-13)  in  Equation  (B-ll)  the 
total  axial  induced  velocity  becomes 

KLIM  JTOT 

W(I,J)  + f—  I I CIRCT(K.MM)  • ,iCTROW(K,tjlM)  (B-U») 

4ttR  k mm 

i B ITOT  JTOT 

+ t—  2 2 £ CIRC(I,M)  ‘ GCB(I.MM) 

4nR  IB  I MM 

Rearranging  Equation  (IB-lU)  to  combine  like  circulations  leads  to 
i B ITOT  JTOT 

W(I,J)  = T—r  E £ I CIRC(lI  ,MM)  * GCDIM(II,MM)  (B-15) 

4ttR  IB  II  MM 

where  GCDBl(lI,MM)  = GCTR0W(II+1,MM)  - GCTRCW(II,MM) 

+ GCT(II,MM) 

Substitution  of  Equation  (B-15)  in  Equation  (B-8)  results  in  a matrix 
equation  in  which  the  only  ’unknowns  are  the  bound  circulation,  CIRC(l,J) 

[CIRC]  = [CONST]  + [GCMAI]  [CIRC]  (B-l6) 


where  GCMAT  = final  matrix  influence  coefficients 

CONST  = matrix  constants  which  are  functions  of  the  following 
known  quantities:  R,  UT,  AT,  C,  and  blade  twist. 


Solving  for  CIRC: 


[CIRC]  = [CONST]  [l-GCMAT] 


(B-17) 
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APPENDIX  C 


GENERALIZED  WAKE  EQUATIONS  FOR  HOVERING  ROTOR 


The  generalized  wake  equations  for  a hovering  rotor,  based  on  the 
model  rotor  experimental  data  of  Reference  8,  are  presented  below.  These 
equations  have  been  taken  directly  from  Reference  8,  where  they  are  de- 
scribed in  greater  detail.  The  flow  visualization  data  in  the  near 
stable  wake  region  of  the  rotor  was  used  in  the  generalization,  and  thus 
the  equations  are  not  intended  to  accurately  represent  the  unstable  far 
wake.  However,  for  performance  calculations  at  the  rotor,  the  wake  in- 
stability generally  occurs  sufficiently  distant  from  the  rotor  to  allow 
the  far  wake  to  be  approximated  by  the  generalized  wake  equations. 

The  generalized  wake  equations  (applicable  in  the  stable  wake  region) 
are  s'mmarized  below. 


Tip  Vortex  Axial  Coordinates: 


k,  \f/w  for  o < <//w  < 


(c-l) 


(ZtU  2tt  + k2  (^w  - 2Jr)  for  - 2J 

Tw  b 

kl  = -0.25  (cT/cr  + 0.001  ff,  J 

k2=-(i.4i  +o.oi4i 

Tip  Vortex  Radial  Coordinates : 

r = A + (I  -A)e~X*w 

A = 0.78  (near  woke) 

X = 0.145  + 27  CT 
Vortex  Sheet  Axial  Coordinates : 


(zr)f  = i 


K|?ii  for  0 - - T 

K-i)¥ 


K..  , = -2.2 

'r*l 


Ct/2 


(C-2) 

(C-3) 

(C-U) 


(c-5) 

(c-6) 

(C-7) 

(c-8) 

(C-9) 

(c-io) 

(c-11) 


BBBGBSB3 


■ 'I111-  .PNWftlSHW-HI 


zr:0 


o for  o < 

(K2f=0  )(*w~ 


7 r 
2 


(C-12) 

(C-13) 


9 , 


K?-  « 
2r:0 


deg 

128 


(° 45  «.d„  + 18 ) I 


(c-D.) 


.‘i:.ce  '.he  cress  sections  of  the  vortex  sheet  from  each  blade  are 
essentially  lines  rather  than  discrete  points,  as  is  the  case  for  the 
tip  vertex  errss  sections,  it  was  assumed  that  the  vortex  sheet  cross 
sections  are  linear,  ar.i  the  axial  position  of  a vortex  sheet  at  a given 
atimuth  can  :e  defined  :y  two  points.  For  simplicity,  the  two  points 
selecte  : were  the  imaginary  extensions  of  the  cross  section  to  r = 0 at 
one  end  ar.i  to  r = 1.0  at  the  other  end  (equations  C-c  throw; h C-lU). 
ihese  two  points  f-sta:lish  the  intercept  at  the  axis  of  rotation  and  the 
slope  of  the  vortex  sheet  at  a particular  wake  azimuth.  The  axial  co- 
ordinatesof  a specific  print  or.  a vortex  filament  forming  part  of  a vor- 
tex sheet  representation  are  obtained  in  the  computer  program  from  this 
intercept,  slope  nr, : the  radial  co<  rdi.nate  of  the  point  defined  below. 

Vortex  Th*  *-t  r.a  iiol  door iinates : 

I4-  was  assumed  that  the  ra  linl  position  of  a point  on  a vortex 
filament  forming  part  of  the  vortex  sheet  re-presentation  (point  A'  is 
linearly  proportional  to  the  radial  coordinate  of  the  intersection  of 
tin-  vortex  sheet  with  the  outer  vortex  sheet  boundary  (point  H)  . With 
the  • xception  of  the  immediate*  vicinity  of  the  blade  (»rw  < 2-rr/V'  , the 
vortex  sheet  boundary  was  assumed  to  be  equivalent  to  the  boundary  formed 
by  the  locus  centers  of  the  tip  vortex  cross  sections  (equivalently,  tip 
vortex  streamline  . For  * less  than  the  blade  spacing,  2-n/b,  the 
boundary  was  faired  from  the  point  of  maximum  circulation  on  the  blade 
(point  Th  to  the  tip  vortex  boundary.  The  constant  of  proportionality 
was  assumed  to  be  the  ratio  of  the  rad.’al  position  of  the  origination  of 
the  vortex  filament  at  the  blade  (point  ~C)  to  the  radial  position  of  the 
vortex  sheet  boundary  at  the  blade  (point  D) . That  is, 


l-'l 


TIP  VORTEX 
CROSS 
SECTION 


J 


5 

(trajectory  i 

OF  TIP 

VORTEX  AND 
OUTER 
VORTEX 
SHEET 
BOUNDARY 


VORTEX 
SHFET 
EXTENSION 
TO  ~ = 1 
(FOR  WAKE 
EQUATIONS) 


FILAMENTS  REPRESENTING  VORTEX  SHEET 


Figure  C-l.  Cross  Section  of  the  Generalized  Wake  Model 


APPENDIX  D 


PRIMARY  EQUATIONS  FOR  THE  UTRC 
PRESCRIBED  WAKE  HOVER  PERFORMANCE  ANALYSIS 


In  order  to  relate  to  the  program  input-output  symbols  and  to 
describe  some  of  the  symbols  generated  within  the  program,  program  symbol 
notation  is  used  in  the  following  equations.  The  equations  for  the  Pre- 
scribed 'Wake  Hover  Performance  Analysis  follow  closely  the  equations  for 
the  Prescribed  Wake  Rotor  Inflow  Analysis  (Appendix  B) , except  that  they 
are  not  time  dependent. 


The  blade  circulation  at  any  blade  segment  I is,  from  the  Kutta- 
Joukowski  relation, 


CIRC(I)  = C(I) *V(I) -CL(l)/2 


(D-l) 


where 


CIRC(l)  = local  circulation,  ft2/sec 

C(I)  = local  blade  chord,  ft 

V(l)  = local  resultant  velocity,  fps 

CL(l)  = local  lift  coefficient 


The  local  resultant  velocity  V(l)  can  be  expressed  as 
V(I)2  = (VPSI(I)  + VTI(I))2  + W(I)2 


(D-2) 


where 

VPSI(I)  = OMGR-RSC(l)  (D-3) 


and  OMGR  = Rotor  tip  speed,  fps 

RSC(l')  = Blade  segment  radial  position 
VT(l)  = Induced  tangential  velocity,  fps 
W(l)  = Induced  axial  velocity,  fps  (also  symbolized  as  VZI  in 
program) . 

For  helicopter  rotors  (versus  propellers),  VT(l)  maybe  neglected.  Also, 
since  W(l)  is  small  in  the  important  blade  tip  region  relative  to  the 
rotational  velocity,  it  is  assumed  that 


V(I)  = OMGR-RSC(I) 


(B-4) 


\ 
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so  that  a linear  closed  form  solution  can  be  obtained.  Also,  in  order  to 
linearize  the  solution,  the  lift  coefficient  is  expressed  as 


I 

l 


CL( I)  = AA(I) -ALPHAE(I) 


(D-5) 


where  AA(l)  = local  lift  curve  slope 

ALPHAE(l)  ■ local  blade  angle  of  attack 

In  stall,  the  maximum  lift  coefficient  is  limited  in  the  program  to 
CLMAX(ll,  which  is  determined  from  input  values  of  ALPHAE  (ALPMAX)  beyond 
which  stall  is  assumed.  The  angle  of  attack  below  stall  is, 

ALPHA  E(I)  = THETA  ( I ) + PHIl(l)  (D-6) 

where  THETA ( 1)  = blade  pitch  angle 
PHII(I''  = blade  inflow  angle 

The  blade  inflow  angle,  PHIl(l)  is 

PHI  I = tan"1  (W(l)/VPSl(l))  (D-7'i 

If  the  tangential  induced  velocity  is  again  neglected,  to  linearize  the 
solution,  it  is  assumed  that  PHII(I)  is  a small  angle  and  thus, 

PHII(I)  = tan  PHII(I)  (D-8) 

Substituting  Equations  (D-E)  through  (D-8)  in  Equation  (D-l)  , the  circu- 
lation equation  becomes 

CIRC( I'  = C(l) -AA(I) • (VPSI(I) -THETA (I)  + W)  (D-9) 

Since  there  are  ITOTI  blade  segments,  there  are  ITOTI  equations  of  the 
form  of  Equation  (D-9^ . This  set  of  ITOTI  equations  can  be  written  in 
matrix  notation  as: 

[CIRC( I) 1 = [C(l)-AA(l) •(  VPS I ( I ) • THETA ( I ) + W(l))]  (D-10) 

The  axial  induced  velocity  in  Equations  (D-9)  and  (D-10)  is  ex- 
pressed, using  the  Biot-Savart  law,  as  a function  of  the  geometry  of  the 
trailing  wake  filaments  and  their  circulation  strengths.  The  axial 
velocity  induced  by  a given  wake  segment  (subscripts  K,  JJJ,  N)  on  trail- 
ing filament  K at  a given  blade  segment  (l,JJ)  is  expressed  in  terms  of 
the  Biot-Savart  law  as 
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W(K,N,JJJ)  = (ClRCT(K)/4nR)-GC  (K,N,JJJ) 


(D-ll) 


CIRCT(K)  = circulation  strength  of  Kth  filament,  ft^/sec 
GC(K = geometric  influence  coefficient  from  the  Biot- 
Savart  law 


The  nondimensional  geometric  coefficient,  GC,  is  a function  of  the  pre- 
scribed wake  cylindrical  coordinates  (RW,  ZW,  and  PHI)  only.  The  lengthy 
expressions  for  GC  may  be  obtained  from  the  program  listing.  Since  the 
blade  circulation  is  assumed  constant  with  time  for  the  hovering  condi- 
tion, each  wake  filament  has  a constant  strength  along  its  length.  The 
total  velocity  due  to  the  Kth  filament  at  the  blade  point  (l,JJ)  is  the 
sum  of  Equation  (D-ll)  evaluated  for  each  segment  in  the  Kth  filament  and 
can  be  expressed  as 


W(K)  = (CIRCT(K) /4nR) *GCK( I, J,K) 


where 


GCK(l,.J,K)  = I EGC(K,N,J.TJ) 
N JJJ 


(D-12) 


(D-13) 


for  a single  K.  Equation  (D-12)  expresses  the  axial  induced  velocity  due 
to  filament  K in  terms  of  the  trailing  filament  circulation.  A summation 
of  the  induced  velocity  contributions  from  all  the  trailing  wake  fila- 
ments, through  a summation  of  the  results  of  Equation  (D-12),  would  lead 
to  an  expression  for  the  total  axial  induced  velocity  at  a blade  segment 
in  terms  of  unknown  trailing  filament  circulations. 

It  will  now  be  shown  that  the  axial  induced  velocity  can  be  expressed 
as  the  product  of  a matrix  of  geometric  coefficients  representing  the  wake 
geometry  and  the  bound  circulations.  Since  all  blades  are  assumed  to  be 
equally  spaced  azimuthally  around  the  rotor  disc,  the  effect  of  the  bound 
circulation  of  all  the  blades  cancels  at  each  blade.  Therefore,  the 
bound  circulation  does  not  directly  contribute  to  the  induced  velocity 
distribution  at  the  blades.  However,  the  circulation  strength  of  each 
trailing  filament  CIRCT(K)  may  be  expressed  as  the  difference  between 
the  bound  circulation  of  the  two  blade  segments  adjacent  to  its  point  of 
origin  at  the  blade.  Thus 


CIRCT(K)  = CIRC( 1=K)  - CIRC(l=K-l) 


(D-14) 


Substituting  Equation  (D-l4)  in  Equation  (D-12)  and  summing  the  effects 
of  all  filaments  from  all  blades,  the  total  axial  induced  velocity  at  a 
point  (l,J)  is 

W(I)  = (l/4nR)S(CIRC(l=K)  - CIRC( I=K-l) ) -GCK(l,  J ,K)  (D-15) 
K 
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Rearranging  Equation  (D-15)  to  combine  like  circulations  leads  to 


b 


W(I)  = 1/4ttRTGCDIM(T,II) -CIRCdl)  (D-l6) 

II 

where  GCDIM(l,II)  = GCK  (l,J,K  = II  +1)  - GCK(l,J,K  = II)  (D-17) 

(summed  on  blade  J values  for  influence  of  wakes 
from  all  blades) 

II  = dummy  subscript  for  I' 

The  azimuthal  subscript,  JJ,  was  eliminated  from  W due  to  the  azimuthal 
independence  of  the  induced  velocity. 

Substitution  of  Equation  (D-l6')  in  matrix  Equation  (D-10)  results  in 
a matrix  equation  in  which  the  only  unknowns  are  the  bound  circulations, 
CTRC(I) : 


[CIRC]  = [CONST]  - tGCMAT]  [CIRC]  (D-l8) 

where  GCMAT  = final  matrix  influence  coefficients  which  are  functions 
of  the  geometric  coefficients  and  input  constants. 

CONST  = matrix  constants  which  are  only  functions  of  the  fol- 
lowing known  quantities:  R,  UT,  AT,  C,  and  blade  twist. 

Solving  for  [CIRC]: 


[CIRC]  = [CONST]  {1  + GCMAT]'1  (D-19) 

A Gauss-Seidel  numerical  procedure  for  evaluating  this  circulation  matrix 
is  included  in  the  computer  program.  Following  the  circulation  solution 
the  induced  velocity  distribution  is  obtained  from  Equation  (D-l6) . 

The  above  equations  are  the  fundamental  equations  in  the  analysis 
for  unstalled,  potential  flow  conditions.  Modifications  providing  for 
blade  stall  and  finite  vortex  core  size  are  included  in  the  computer 
program. 
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LIST  OF  SYMBOT.S 


blade  coning  or  precone  angle,  deg 


longitldinal  flapping  angle  for  the  first  harmonic; 

of  ^r  'T*  te™  “ F0Urier  Series  mansion 

— to  the 


wing  area,  ft2 


blade  aspect  ratio,  AR  R/c 


cosine  component  of  first  hannonic  cyclic  pitch; 

~ - =s  - 


number  of  blades 


ojn-sI„Ttem  in  FirSt  CyCliC  pit^  -efficient 

sin*  term  in  Fourier  series  expansion  of  blade  nitrh 

angle  with  respect  to  the  plane  nonsal  to  the  rotor  shaft, 


blade  chord,  in. 

rotor  lift  coefficient,  L/prrR2(nR)? 
wing  lift  coefficient,  L^/ !pV2  Aw 
rotor  thrust  coefficient,  T/prrR2 ( fJR ) 2 
rotor  lift,  lb 
wing  lift,  lb 

dynamic  pressure,  IpV2,  pSf 


radial  coordinate  from 


z or  z’  axis,  ft 


rotor  radius,  ft 


1 


rotor  thrust,  lb 

momentum  theory  value  of  rotor  induced  velocity, 
absolute  value,  fps  


4F 


Hover:  vD  = CB 


Forward  Flight:  vc 


HR  Cf 


vx>  vy  > v z 


V’  v>  vz’ 


VX>  Vy,  Vz 


VX’>  Vy ' , V, 


local  instantaneous  induced  velocity  component  in  the 

x,  y,  or  z directions,  positive  in  positive  axis 
direction,  fps 

local , instantaneous  induced  velocity  component  in  the 

x 5 y , or  z'  directions,  positive  in  positive  axis 
direction,  fps 

local  time-averaged  induced  velocity  component  in  the 

x,  y,  or  2 directions,  positive  in  positive  axis  direc- 
tion, fps 

local, time-averaged  induced  velocity  component  in  the 

x 5 y , or  z’  directions,  positive  in  positive  axis 
direction,  fps 

velocity  component  induced  at  rotor  by  wing,  normal  to 
rotor  disk,  positive  up,  fps 


free-stream  velocity,  fps 


longitudinal  coordinate  from  rotor  hub  center  in  plane 
parallel  to  rotor  disk,  positive  in  downstream  direction, 

roto^h”?  rrdlnate  “ "in<:  “U  “^sured  fom 

center,  positive  in  downstream  direction,  ft 

tfroTi  C°0r?lnate.fro”  rotor  hub  center  in  plane  parallel 
isk,  positive  toward  advancing  side  of  rotor,  ft 

lateral  coordinate  from  rotor  hub  center  in  wind  axis  sys- 
tern,  positive  toward  advancing  side  of  rotor,  ft 
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coordinate  from  rotor  hub  center  in  direction  normal 
to  rotor  disk,  positive  up,  ft 

z/rotor  radius,  R 

vertical  coordinate  from  rotor  hub  center  normal  to 
wind  axis,  positive  up,  ft 

coordinate  from  center  of,  and  normal  to,  the  rotor 
tip- path  plane,  positive  up,  ft 

vertical  coordinate  from  local  point  on  a blade,  positive 
up,  ft 

same  as  z except  measured  from  tip  of  blade 

rotor  tip-path-plane  angle  of  attack.  Angle  between 
free- stream  velocity  and  tip-path-plane 

local  blade  circulation  or  circulation  of  a wake  vortex 
element,  ft 2/ sec 

spanwise  average  of  downwash  angles,  rad,  e sb  vz/Vcos a 

theoretical  downwash  angle  at  center  of  rotor  as 
defined  in  Reference  27,  e0  — v0/Vcosor 

blade  linear  twist  rate  from  center  of  rotation  to  tip, 
deg  (normally  negative) 

blade  collective  pitch  angle  measured  at  r = 0.75  R} 
deg 


rotor  inflow  ratio,  X = (Vsinor  - v0)/flR 

rotor  advance  ratio,  ^ = Vcoso  /GR  (or  V/GR  as  indicated) 

air  density,  Ib-sec^/ft*4 


rotor  solidity  ratio, 


CT 


bc/rrR 
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